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Abstract 
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dimensional field theory are discussed, among them the running of gauge couplings, the 
spectrum of mesons in Higgs and Coulomb branches, entanglement entropy and the non- 
relativistic version of our string backgrounds. 
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1 Introduction 



Since its inception p], the large N c expansion has given many valuable insights into QCD 
dynamics beyond perturbation theory. Indeed, long ago 't Hooft showed, by using this 
expansion, that two-dimensional models are very good laboratories in which one can learn 
about various aspects of four-dimensional field theories, like the spectrum of bound states 
[2]. Besides, the large N c expansion sheds light on other aspects, like the connection between 
the Skyrme model and QCD and the qualitative understanding of the Zweig rule [3]. 

The number of colors is not the only available parameter in QCD-like theories. In the 
real world, and specially at energies above the masses of the heavy quarks, the ratio Nf/N c 
between the number of flavors and colors becomes of order one and, then, it seems natural 
to investigate the way in which a large number of flavors modifies the dynamics of the large 
N c field theory. From a field theory point of view, these ideas were introduced in [4]. 

In this paper we consider these kind of problems from the perspective given by the 
AdS / CFT correspondence [5] , its refinements [6] , [7] and extensions [H [9] . In this approach 
the addition of matter degrees of freedom can be performed by including extra (flavor) branes 
in the description [TO] and one has to deal with a system of gravity coupled to brane sources 
|f lj . When the number of flavors Nf is small compared with the number of colors N c one 
can treat the flavor branes as probes in the geometry originated by the color branes. This is 
the so-called quenched approximation. However, if Nf ~ N c the backreaction of the flavor 
branes on the geometry cannot be neglected anymore and one must obtain a new background 
that incorporates the effect of the fundamental matter. 

We will apply these ideas in order to study the case of supersymmetric field theories 
in two spacetime dimensions. This paper represents a first step in this direction as we 
will concentrate in the case of field theories with M = (4, 4) SUSY. The analysis of two- 
dimensional theories with M = (2,2) and M = (1, 1) supersymmetry is in preparation and 
will be reported elsewhere. 

We will proceed by constructing string duals to Af = (4, 4) gauge theories. Our back- 
grounds will be based on D3-branes wrapping two-cycles of a Calabi-Yau manifold of com- 
plex dimension two (leading to a 2d field theory at low energies) or on D1/D5 systems of 
fractional branes. We will explore our backgrounds with "probe branes" to learn about var- 
ious aspects of the dual QFT, most notably, the running of couplings and the "quenched" 
spectrum of mesons (in the Higgs or Coulomb branches). Then, following the line of research 
presented in the papers [12]- [21], we will add a (large) number of flavor branes to the ge- 
ometries mentioned above, which correspond to flavor degrees of freedom in the dual QFT. 
We will then find string backgrounds dual to M = (4, 4) QFT's in two dimensions with large 
N c ,Nf. Again, we will present some checks of the field theory-string theory matching and 
obtain some non-perturbative predictions from our backgrounds. 

In the case of backgrounds generated by D3-branes wrapped on a two-cycle, we will adopt 
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an ansatz for the metric and five-form which contains functions depending on two radial 
coordinates. By imposing the preservation of eight supersymmetries we will find a system of 
partial differential equations. These equations look rather simple, but are difficult to solve 
in general. However, we will be able to find some particular analytic solutions. It turns out 
that the most interesting of these solutions is provided by the analysis of five-dimensional 
gauged supergravity, where supersymmetry is realized by means of the so-called topological 
twist [22] • Indeed, we will verify that this solution encodes non-trivial information on the 
d — 2, J\f = (4, 4) gauge theory. Moreover, we will be able to modify the BPS equations 
to include the backreaction of the flavor branes, which in this case are D3-branes extended 
along the non-compact directions of the normal bundle of the cycle where the color branes 
are wrapped. The corresponding solutions, which in this case are numerical, also pass several 
non-trivial tests. 

An alternative string description of the d — 2, jV = (4, 4) gauge theories is obtained by 
considering a system of fractional Dl-branes on the orbifold R 1 ' 5 x C 2 /Z 2 . In this case one 
proceeds by solving directly the second-order equations of motion for the fields of type IIB 
supergravity, including the ones originating from the twisted sector of the orbifold theory. 
In this approach the flavor branes are D5-branes that wrap completely the C 2 /Z 2 orbifold. 
We will verify that this solution also matches the field theory results. 

The organization of the rest of this paper is the following. In section [2] we will introduce 
our setup of D3-branes wrapped on a two-cycle. The detailed derivation of the system of 
BPS equations are presented in appendix [A] The gauged supergravity solutions of this BPS 
system are obtained in appendix [B] and analyzed at the end of section [21 In section [3] we 
study the addition of flavor to our solution and we obtain the corresponding backreacted 
geometry. In section H] we first review the field content of the d — 2, J\f — (4, 4) gauge theories 
and the calculation of their one-loop beta functions. By performing a probe calculation we 
check that our backgrounds match these field theory results. We also analyze in this section 
the Higgs branch of the theory, based on the study of the supersymmetric embeddings of 
D3-branes performed in appendix We finish section H] by analyzing the entanglement 
entropy of our model in the UV region, following the recent proposal of refs. [23] and [24J. 

Section [5] deals with the study of the mass spectra of mesons for our theory, starting with 
the case of the Coulomb branch in the quenched approximation. We also discuss the effect of 
the backreaction, as well as the spectra for the Higgs branch. In section [6] we study additional 
solutions of the unflavored BPS system. In particular we show that our equations also admit 
a solution in which the geometry has an AdS% factor. We also obtain the background dual 
to a non-relativistic system, which is generated from our solution by means of a pair of 
T-dualities, combined with a suitable coordinate shift. Section [7] is devoted to the analysis 
of the system of D1/D5 fractional branes. Finally, in section [8] we summarize our results and 
discuss some of their possible extensions. 
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2 The supergravity dual of the (4,4) gauge theories 



Let us consider the background of type IIB supergravity created by a stack of N c D3-branes 
wrapped on a two-cycle C 2 of a Calabi-Yau (CY) cone of complex dimension two according 
to the brane setup: 



CY 2 





R 1,1 


S 2 


N 2 


R 4 


D3 

























where S 2 represent the directions of the two-cycle (which is topologically a two-sphere) and 
N2 are the directions of the normal bundle to C2 ■ Notice also that the symbols "— " and "■" 
represent respectively unwrapped worldvolume directions and transverse directions, while a 
circle denotes wrapped directions. 

Let us parameterize C2 by means of two angular coordinates (9, 0) and let a be the radial 
coordinate of the CY cone. Notice that in this setup there is another radial coordinate 
p, which represents the distance along R 4 , the directions orthogonal to both the D3-brane 
worldvolume and the CY cone. The ansatz for the string frame metric which we will adopt 
is the following: 

iH \dx\ l + A; ( d9 2 + sm 2 6d(t) 2 )] + 



ds 



St 



- {da 2 + a 2 

z 



dip + cos 



+ dp 2 + p 2 dQA 



where m is a constant with units of mass which, for convenience, we will take as: 



— = y/4ng 8 N c 



a 



(2.1) 



(2.2) 



with g s and a' being respectively the string coupling constant and the Regge slope of su- 
perstring theory. In eq. ( 12. 11) dxf^ denotes the two-dimensional Minkowski metric for the 
coordinates x°, x 1 and the range of the variables 9, <fr, ip, p and a is the following: 



< 9 < 7T 



< 0, ip <2n , < p , a < 00 . 



(2.3) 



Moreover, the function z (which controls the size of the cycle) and the warp factor H should 
be considered as functions of the two holographic variables (p, a): 



H = H(p,a) 



z(p,a) 



(2.4) 



while if}, which is a coordinate of the CY, is fibered over the cycle Moreover, dQ 2 is 
the metric of a unit three-sphere, which we will parameterize in terms of three angles «j 
(i = 1, 2, 3) as: 

dfi 2 , = da\ + sin 2 a,\ ( da 2 , + sin 2 02^0:3) . (2.5) 
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As in any background created by D3-branes, our solution should be endowed with a a self- 
dual RR five-form F$, that we write as: 

F 5 = T h + *F 5 , (2.6) 

where JF 5 can be represented in terms of a four- form potential C4 as JF 5 = <iC 4 . We shall 
adopt the following ansatz for C 4 : 

C4 = g cu 3 A (dip + cos 9 d(j)) , g = g(p,a), (2.7) 

where g(p, a) is a new function and u 3 is the volume element of the three-sphere: 

to 3 = sin 2 cti sin a 2 da,\ A A da 3 . (2.8) 

To determine the functions H, z and g entering our ansatz we will impose that our back- 
ground preserves eight supersymmetries, which is the right number of SUSYs that the super- 
gravity dual of a d = 2, Af = (4, 4) gauge theory must leave unbroken. The detailed analysis 
is performed in appendix [A] It turns out that, in order to preserve this required amount of 
supersymmetry, the three functions H, z and g must satisfy the following system of partial 
differential equations (PDEs): 

2 3 / 

m g = p z , 

2 rr Zi 

m H = — , 
a 

g' = -ap 3 H , 

9= ^H'-^m"gH, (2.9) 

z z l 

where we have denoted: 

' = d p , -=d a . (2.10) 

Actually, the four BPS equations (12.91) are not all independent. Indeed, one can check 
that, for example, the last equation in the system (12. 9p is a consequence of the others. 
Furthermore, one can combine the first-order BPS equations (12. 9p to obtain the following 
second-order PDE equation for z(p,a): 

pz(z — a z) = cr(pz 2 + pz" + 3z' ) . (2-H) 

Moreover, once the function z(p, a) is known, one can obtain, by using the first two equations 
in (12.91) . the values of the other two functions of our ansatz, namely H(p, a) and g(p, cr). 
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It is easy to verify that, when the system (12.91) holds, the field strength F 5 satisfies the 
Bianchi identity dF$ = 0. Indeed, one can check this fact directly by using the equations in 
(12. 9p (see appendix [AJ or, alternatively, one can verify that F$ can be represented in terms 
of a four-form potential C4 as F 5 = dC^. The actual expression of C4 can be taken as: 

C 4 = jw 3 A (dif) + cos 6d<p) + dx° A dx 1 A 

where uo 2 is the volume element of the (8, (/>) two-sphere, namely: 

lu 2 = sin BdB A d<f> . (2.13) 

One can also verify that the system (12.91) ensures the fulfillment of the second-order Einstein 
equations. The details of this verification are given in appendix [A] 



a 



. m, 2 H 



lu 2 da A (dip + cos 6d(p) , (2.12) 



2.1 Integration of the BPS system 

The system (12. 9ft can be easily integrated for g and z when a = and p varies. Indeed, it 
follows from the last two equations in (12. 9p that the five-form function g is independent of p 
when o = 0. Let us write: 

g(p,a = 0)=g , (2.14) 

with go being a constant which should be determined by flux quantization of the RR five-form 
(see below). Then, the first equation in (I2.9P for a = can be written as: 

^'(P,0) = ^, (2.15) 
P 3 



which can be integrated as: 



m g 

z(p,0) = 5 h constant . (2.16) 

2p z 



In order to get an explicit solution of the BPS equation (I2.9p for all values of the (a, p) 
coordinates, it is interesting to notice that the brane setup analyzed here can be studied in the 
framework of five- dimensional gauged supergravity. Actually, the corresponding topological 
twisting was studied in [22] and is summarized in appendix [Bl One of the advantages of 
using gauged supergravity is the fact that the non-trivial fibering of the ip coordinate in the 
metric comes up naturally in the uplifting from five to ten dimensions as a consequence of 
the twisting. Also, the RR five-form is given by the uplifting formulae of ref. [25]. The 
ansatz of gauged supergravity contains three functions which depend on a single holographic 
coordinate (see appendix [B]). The corresponding BPS equations are a system of ordinary 
differential equations which, as shown in appendix [Bl can be analytically integrated. After 
performing a suitable change of variables the uplifted metric and five-form can be written as 
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in (12. ip and (I2.6l) - (l2.7p . This is worked out in detail in appendix [B] Let us summarize here 
the main results. With this purpose, we define the function T(z) as follows: 

r(z) = z* + (z* - z) [ \og{z* - z) + k] , (2.17) 

where and k are constant^. The function T(z) defined in (12.171) can be used to determine 
implicitly z as a function of (p, a). Indeed, it is shown in appendix |B1 that z(p, a) satisfies: 

(2.18) 



P 



r z) 



z*- z 



2m? 



where m is the same constant as in (12.21) . Notice that, by taking a = in (12.181) one can 
immediately solve for z, with the result: 



z(p,a = 0) 



1 



2m 2 p 2 

This result is consistent with (12.161) if the constant go is taken to be: 

1 



9o 



rrr 



Ang s (a') 2 N c . 



(2.19) 



(2.20) 



One can check that, non-trivially, the function z(p, a) defined implicitly in (12.181) satisfies 
the second-order PDE (12. lip (see appendix IB]) . 

From the gauged supergravity approach of appendix [B], one can also find the function 
g(p, a) in terms of z(p, a). One obtains: 



9 = —o 



^z* - z)p* 



m 2 p 2 + 



(2.21) 



As a check of (12.211) it is easy to verify using (12. 191) that g(p,cr = 0) = g = m 4 . Moreover, 
one can also obtain the expression for the warp factor H, namely: 



H 



2z(z* - z) 



p 2 + 



a" 



T(z) 



(2.22) 



In particular we can now obtain the value of the warp factor H at o = 0. One obtains: 

1 2m 2 p 2 z* — 1 



#0,(7 = 0) 



m 4 p 4 2m 2 p 2 z* — log(2m 2 p 2 ) + k 



(2.23) 



Let us now study the solution for arbitrary values of p and a. First of all, it is clear 
from (12.181) that, when either p or a are large, the function z approaches the constant value 



In all the numerical calculations performed in this paper we take k = . 
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Figure 1: From left to right and from top to bottom plots of z, H and gin 4 versus the 
dimensionless variables pm and am. We are setting = 5 and k = 0. 

z = z*. Moreover, by solving numerically the implicit relation ( 12.181) one can obtain the 
function z(p, a) and, then, by using this result in (12.221) and (I2.2ip one can obtain H and 
g. The result of this numerical analysis has been plotted in figure [TJ One important thing 
that can be observed in this result is that the function z becomes negative when p and a are 
small enough and, therefore, the supergravity solution is inconsistent in this region. This 
phenomenon is related to the so-called enhangon mechanism [26J. Indeed, we will explicitly 
show in subsection 14.11 that when z vanishes the color brane probes become tensionless, 
signalling the appearance of new massless degrees of freedom in the IR. One can estimate 
the scale at which this mechanism starts by computing the value of p at which z vanishes 
for a = 0. From ( 12.191) we obtain that this value of p is given by: 



Thus, the asymptotic value of z also controls the IR scale of our solution. 

In spite of the fact that we cannot solve analytically for z(p, a) in the implicit equation 
( I2.18p . one can solve this equation by means of an expansion around the constant asymptotic 
value z = z*. Notice that, at zeroth order in this expansion, T(z*) rs z*. Using this result to 
evaluate the left-hand side of (12 . 181) one finds that, keeping the first non-trivial term, z(p, a) 



1 



(2.24) 
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can be approximated to be: 



z(p, a) 



2m 2 a 2 + z* p 2 



(2.25) 



Notice that this expression gives the exact result for a = (see eq. (12.191) ). Moreover, it gives 
a rather good approximation for o > 0. Actually, when one compares the function z(p, a) 
obtained by a numerical solution of (12.181) with the one written in (12.251) . one realizes that 
the two functions only differ significantly when we approach the enhangon point. Moreover, 
by substituting z(p, a) as given by (12.251) in the right-hand side of eqs. (12.211) and ( 12.221) 
one gets a very good approximation to the functions g and H. These expressions greatly 
simplify if we assume that p or a are large and we keep the leading terms. In this case, one 
gets: 



rrr 



a 2 + z*p^ 



2 ' 



P' 



rrr 



a 2 + z*p 2 



2 ' 



(2.26) 



which are good estimates far from the enhangon point. This analysis suggests that in the 
UV region the particular combination a 2 + p 2 of the radial variables plays a relevant role. 
Accordingly, we define new variables u and a as follows: 



u 



o 2 + z*p 2 



tana 



a 



z*p 



7T 

< a < - . 
~ ~ 2 



(2.27) 



The functions H and z of the metric will generally depend on both coordinates u and a: 

H = H(u,a) , z = z(u,a) . (2.28) 

The metric in the new coordinates takes the form: 



ds 2 



H-2 



dx\ 1 H {d0 2 + sin 2 9 1 



sin 2 a cos 2 a 



+ 



du 2 + 



2 / cos 2 a sin 2 a \ 2 A A / 1 1 

+ u 1 da + 2u sin a cos a 



z 



z 



duda + 



u 



u 



H sin a (rf^ + cos ^rf^) + — cos a dQ 



z* 



(2.29) 



It follows from (I2.25|) and (12.261) that, when u is large, the functions H and z become 
independent of a and are given by: 



z ^ z± 



H 



m 4 u 4 



(2.30) 
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Using these values for z and H in the metric ansatz , we get the following: 



ds 



uv 



m 



u 



dx\i + ~^(d9 2 + sin 2 9, 



+ 



1 du 2 



+ 



+ — 

mr 



da 2 + sin 2 a (dip + cos 



2 + cos 2 a dfl 2 - 



(2.31) 



The first line in (12.311) is the metric of an AdS§ space, in which two of its directions are 
compactified on an S 2 . The second line of (I2.3ip is the line element of an S 5 fibered over 
the S 2 . Notice that this result is in agreement with the origin of the solution as wrapped 
D3-branes in gauged supergravity, with the fibering being given by the gauge field in this 
approach. Actually, the variables u and a. can be easily identified with the ones used in the 
gauged sugra approach. In this case the UV region corresponds to the large r region in eq. 
(IB. 191) of appendix [Bl From the definition of p and a in (IB.41I) one has the following relation, 
in the UV, between the coordinates (u, a) and the five- dimensional coordinates (r, #): 



u 



a 



a 







rn 



(2.32) 



where a is the integration constant appearing in flB.19j) . Notice that, with this identification, 
the metric flB.28j) becomes for large r exactly the one written in (I2.3ip . 



3 The dual of the (4,4) theory with flavor 

In the previous section we have succeeded in finding the gravity dual of M = (4, 4) super 
Yang-Mills theory in two space-time dimensions with a vector multiplet in the adjoint repre- 
sentation of the gauge group. The main purpose of the present section is the extension of this 
result to include matter hypermultiplets in the fundamental representation. We will follow 
the now standard procedure of adding flavor branes to the setup of color branes discussed 
in section [2j These flavor branes add a new open string sector to the theory, which can 
be interpreted as the gravity dual of the addition of matter hypermultiplets on the gauge 
theory side. The flavor branes should fill the Minkowski spacetime and should extend along 
some holographic direction while wrapping some internal submanifold of the Calabi-Yau 
manifold. We will show below that the appropriate flavor branes for our case are D3-branes 
extended along the Minkowski space-time directions well as the directions (a, ip) 

of the normal bundle of the Calabi-Yau, according to the array: 



CY 2 





R 1,1 


S 2 


N 2 


R 4 


N c D3 (color) 
























N f D3 (flavor) 
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The Nf D3-branes in this array are located at a fixed point on the S 2 sphere parameterized 
by the angles (8, 0) and at a fixed point in the transverse R 4 . Actually, the value pQ 
of the coordinate p of the flavor branes represents the distance between the two sets of 
branes and, according to the radius-energy relation ( 14.211) . is related to the mass itlq of the 
hypermultiplets as follows: 

m Q = — ■ (3.1) 

One of the arguments in favor of considering these configurations as the right ones to add 
flavor to the Af = (4,4) gauge theory is the fact that this setup preserves the same supersym- 
metries as that considered in section [2j This statement is verified in appendix [C] by using 
the kappa symmetry of the DBI action. Notice also that, as expected on general grounds, 
the flavor D3-branes of our setup are extended along the non-compact direction of the CY 
cone. In general, if X M denote ten-dimensional coordinates, the D3-brane embedding will 
be characterized by a set of functions X M (£ a ), where £ a (a = 0, ■ ■ • , 3) is a system of world- 
volume coordinates of the D3-brane. In what follows we shall choose the following set of 
worldvolume coordinates: 

C = {x\x\^,a) , (3.2) 

and our embedding will be characterized by having a constant value of the remaining ten- 
dimensional coordinates, namely: 

p = pq , 9 , , «j = constant . (3.3) 

Notice the particular ordering we are using for the coordinates £ a in ( 13. 21) . This ordering 
determines the orientation of the worldvolume of the flavor branes, which in our case differs 
from the one induced by the ten-dimensional background. 

When Nf << N c one can neglect the backreaction of the flavor branes on the geometry 
and, thus, one can treat them as probes. This is the so-called quenched approximation, which 
corresponds, in the field theory side, to suppressing quark loops (that contribute with powers 
of 1/N C ) in the 't Hooft large N c expansion. By analyzing the normalizable fluctuations of 
these brane probes one can extract the meson spectrum of the model (see [11] for a review). 
We will postpone this analysis until subsection 15.11 In the remainder of this section we will 
study how to incorporate the effect of the backreaction in our brane setup. 



3.1 Including the backreaction 

When the number of flavors Nf is of the same order as the number of colors N c the backre- 
action of the flavor branes cannot be neglected anymore and one is led to consider the full 
coupled gravity plus branes system. Including the backreaction is the analogue, on the field 
theory side, of considering the effects of quark loops, which are suppressed in the 't Hooft 

2 This does not happen if one adds D7-branes as flavor branes. 
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large N c limit. In what follows in this subsection we will assume that Nf is large and Nf/N c 
is fixed. To construct the corresponding supergravity dual we will follow the approach pio- 
neered in ref. [12] and we shall consider a suitable continuous distribution of flavor branes 
(see [27] for a similar analysis in the context of non-critical string theory). This approach 
has been successfully applied to study backgrounds with different amount of supersymmetry 
which are dual to gauge theories in several space-time dimensions [13]- [21]. In section [8] we 
will comment on an interesting subtlety about the smearing and the validity of the DBI+WZ 
action 

The starting point in our analysis is the observation that the embeddings considered at 
the beginning of this section are mutually supersymmetric for any value of the coordinates 
transverse to the flavor D3-brane. Therefore, when Nf — > oo, we can homogeneously dis- 
tribute the Nf flavor branes in some of their transverse directions. Actually, we shall locate 
them at a particular value p = pq of the p coordinate, and we will smear them along the 
angular directions (9, 4>) as well as along the three-angles of the three-sphere. 

The action of a stack of Nf D3-branes is given by: 

S flavor = Sdbi + Sy/z j (3.4) 
where the DBI and WZ terms are given by: 

Sdbi = -T 3 ^ / yf - det G 4 , S wz = T 3 ^ / C 4 , (3.5) 

N f Nt 

with G4 being the induced metric on the worldvolume M.± and the hat over the RR four-form 
potential denotes its pullback to M.±. The smearing procedure^ amounts to performing the 
following substitution in the WZ term in (13.51) : 



N f 



C 4 -> / OAC 4 , (3.6) 

Ma J M10 



where Q is a six-form proportional to the volume form of the space transverse to the world- 
volume of the flavor brane: 



N 



« = Vol(;y 6 ) , (3.7) 



57T° 

with Vol (y e ) being: 

Vol (y 6 ) = 5{p - p Q ) dpAu 3 Auj 2 . (3. 



3 Notice that we add brane sources via the DBI+WZ action. These are the branes we smear. This is 
different from what is studied in |28j and hence the comments in that paper do not apply to our setups. On 
the contrary, our approach is closer to [29], where a general analysis of the eqs. of motion was performed 
and shown to be always consistent. We are grateful to Ingo Kirsch for discussions on this point. 
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The normalization constant in f2 has been chosen to satisfy the normalization condition: 

tt = -N f , (3.9) 



where the minus sign is due to the orientation of the worldvolume. Therefore, after the 
smearing, the WZ term of the action of the flavor branes takes the form: 



S wz = T 3 fiAC 4 . (3.10) 

JMio 

In the DBI part of the action the smearing is performed by means of the substitution: 



J2 [ d 4 i \j- det G A -> I d 10 x V-detG\Q\, (3.11) 



N f 

where I Q is the modulus of the form f2, defined as: 



n 



\ 



1 6 

-fl Ml ... M(i Q Nl ... N6 J~J G MkNk . (3.12) 



6! 

fc=i 



The smeared DBI+WZ action of the flavor branes is thus: 

*S flavor 



-T 3 I d 10 x V- det G | tt | + T 3 f OAC 4 . (3.13) 

J Mio J Mio 

We will assume that the metric of the backreacted background will be still given by the 
ansatz (12.11) . i.e. it can be written in terms of the warp factor H(p,a) and the function 
z(p,a). On the contrary, it is evident from the form of the WZ action in (13.101) that the 
flavor branes act as sources of the RR form F 5 and, as a consequence, they will induce 
a violation of its Bianchi identity in the backreacted geometry. Therefore, to obtain the 
backreacted background we are looking for, we will have to modify our ansatz (I2.6p -( 12T7I) for 
F 5 . It is clear from (I3.10P that the modified Bianchi identity in this case is: 

dF 5 = -2tt g s (a') 2 N f 5(p - p Q ) dp A u 3 A u 2 . (3.14) 

Accordingly, let us represent F 5 as in (12. 6p with being given by: 

F 5 = h- 27Tg s (a , ) 2 N f Q(p~p Q )uj 3 Aiu 2 , (3.15) 

with / 5 such that df$ = 0. We shall represent /s in terms of a potential C4 as /s = dC A , 
where C 4 is parameterized by the function g(p, a) as in (12.71) . 
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Proceeding as in the unflavored case and substituting our ansatz for F 5 in the equations 
for the SUSY variations of the dilatino and gravitino (see appendix EJ), we get the following 
set of BPS equations: 



m 



2 [g - 2ng s (a') 2 N f e(p-p Q )} = p 3 z' 



2 rr ZZ 

m H = — , 
a 



-ap 3 H 



g = ^H' - - 2 Hm 2 [g - 2n g s (a 1 ) 2 N f e(p-p Q )] . (3.16) 

Clearly, when Nf = or p < pq the system (13.161) reduces to (12. 9p . Moreover, one can prove 
that z(p, a) satisfies the following PDE: 

pz(z - az) = a(pz 2 + pz" + 3z' ) + 5(p - p Q ) . (3.17) 

2J\I C m z p z 

Notice that the flavors contribute in (13.171) as a source localized at p = pq. 

The set of projections to be imposed on the Killing spinors in order to arrive at the 
system (I3.16P is just the one written in (IA.7I) . i.e. they are same as in the unflavored case. 
Therefore, any solution of (13.161) preserves eight supersymmetries. Moreover, as in (12.91) . the 
last equation in (13.161) can be derived from the first two. On the other hand it can be checked 
that the equation (13.141) for F 5 is just a consequence of the system (13.161) (see appendix 
We also verify in appendix |X] that the Einstein equations, including the contribution of the 
DBI term of the flavor branes action (I3.13p . are satisfied by any solution of (13.161) . 

As in the unflavored case, we can integrate the function z for a = 0. Indeed, it follows 
from (13.161) that g is independent of p when the variable a vanishes. Let go denote this 
constant value of g. Then for p > pq, we have: 

z'(p,0) = m 2 [g - 2ng s (a') 2 N f j \ , (p > p Q ) . (3.18) 

Using for g the same value as in the unflavored case, the above equation becomes: 



N 



z'(p,0) = 4nm 2 g s (a') 2 N c - ^, (p > p Q ) . (3.19) 



2 



1 



P 3 



When p < pq one simply puts Nf = on the right-hand side of (13.191) and (12.151) is recovered. 
Notice that z'(p, 0) jumps when p passes through the point p = pq. Actually, the slope of 
the z(p, 0) curve becomes negative when Nf > 2N C . We will argue below that this is the 
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Figure 2: Plots of z(p, a) for the backreacted case. On the left we present the solution with 
Nf = N c and on the right that corresponding to Nf = 7N C . On the right graph we have 
zoomed in on the region of interest so that one can easily see that the slope becomes negative 
at p = pq. We have set z* = 5 and pq = 1.5. 



geometric counterpart of the behavior of the beta function in the field theory dual. Let 
us now integrate (13.191) and impose the condition that the solution z(p, 0) is continuous at 
p = Pq (although its first p derivative is not). One gets: 

irm 2 q s (a') 2 , 2ir m 2 a* (a') 2 r N f , i 

z(p,0) = z* y^U- Nf e(p-n) N c -^Q(p- PQ ) ,(3.20) 



2 



where the constant is the same as in ( 12.19p . In order to evaluate z for arbitrary values of p 
and a we must numerically integrate the flavored system (I3.16P . The result of this numerical 
integration in two particular cases is presented in figure EJ To obtain these results we have 
assumed that the solution reduces to the unflavored one for p < P q. We also assumed that 
g(p,cr) is continuous at p = P q and thus (see the first equation in (13.161) ) z'(p,a) has a 
discontinuity at p = pq which is independent of a and given by: 

1 N f 

z( PQ + e,a)- z( PQ ~^°) = -^ qWc - (3-21) 

In the plots of figure [2] it is quite evident the qualitative difference between the cases with 
Nf < 2N C and Nf > 2N C . Indeed, when Nf < 2N C , the function z continues to grow for 
p > P q as we move away from the origin in the (p, a) plane, whereas for Nf > 2N C there 
is a change of behavior at p = pg, where z starts to decrease and, if pg is not too large^, z 
becomes negative and the space ends. 



One can estimate the maximal value of pq such that z becomes negative when Nf > 2N C from the exact 
expression p.20[) of z at a = 0. Indeed, a simple analysis of (|3.20p shows that this maximal value of pq is 

p* , where p* is the IR scale defined in (|2.24l) • 
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4 Matching the field theory 



The two-dimensional gauge theories with Af = (4, 4) supersymmetry can be obtained by 
dimensional reduction of four- dimensional Af = 2 gauge theories 0. As in four-dimensional 
Af = 2 theories, the Af = (4, 4) theories in d = 2 have two massless representations, namely 
the vector multiplet and the hypermultiplet. 

The field content of the d = 2, Af = (4, 4) vector multiplet can be obtained by dimensional 
reduction of the corresponding vector multiplet of the d = 4, Af = 2 theory Recall that the 
latter is composed of a vector field A^, a Dirac spinor A and a complex scalar <p. To obtain 
the dimensional reduction of the vector field one simply distinguishes the components 
with fi = 0, 1 from those with \i = 2, 3. This results in a two-dimensional vector field and two 
real scalars. Moreover, the dimensional reduction of the d = 4 Dirac spinor A gives rise to 
two complex fermions in d = 2 whereas, by dropping the dependence on the coordinates x 2 
and x 3 of the d — A complex scalar <ft, one gets two more real scalars in the two-dimensional 
theory. Therefore, the total content of the d = 2, Af = (4,4) vector multiplet is one vector 
field, two Dirac fermions and four real scalars (which parameterize the Coulomb branch of 
the theory). 

The d = 2, Af = (4, 4) hypermultiplet is similarly obtained from the d = 4, Af — 2 
hypermultiplet. The latter contains a Dirac fermion and two complex scalars which, upon 
dimensional reduction, give rise to two Dirac fermions and four real scalars. These scalars 
parameterize the Higgs branch of the theory 0. 

The R-symmetry of the (4, 4) theory is SU (2) rX SpiniAf); the first part is inherited from the 
R-symmetry of the six dimensional theory with eight supercharges and the Spin(A) ~ SO (4) 
is coming from the dimensional reduction on four directions of the initial six. Notice that 
our background (12. ip realizes the SO (4) explicitly as rotations in the three-sphere whose line 
element is dVL\, while it only realizes explicitly a U{1) subgroup of the SU(2) R , represented 
by traslations in the angle if. 

Given the field theory content described above, it is possible to find the running coupling 
constant in perturbation theory for the Af = (4, 4) super Yang- Mills theory [30]. Here we will 
follow the presentation of the appendix C of [31] , for the particular case of two-dimensional 
theories. The one loop running coupling constant is given by: 



5 Alternatively, they can also be obtained by dimensional reduction of Af — (1, 0) supersymmetric theories 
in six dimensions. 

6 The vector and hyper multiplets of the d = 2, Af = (4, 4) theories can also be regarded as being composed 
of supermultiplets of d = 2, Af = (2, 2) supersymmetry. In terms of the Af — (2, 2) superalgebra the vector 
multiplet contains a chiral multiplet and a twisted chiral multiplet, while the hypermultiplet decomposes 
into two chiral multiplets. 




(4.1) 
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where // is the energy scale and b is a constant which depends on the field content of the 
theory. Let us consider a theory with n s scalar fields, n v vector fields and n f Dirac fermions. 
In this case the value of b is the following: 

n s 2 
b = y c s - 4n v c v + -rifCf , (4.2) 

where cj for I = s,v,f is the normalization constant of the generators of the gauge group 
in the representation of the different fields (tr(T a T b ) = ci5 ab ). With the conventions chosen 
in [3T] one has for the gauge group SU(N C ): 

^fundamental ^_ ^adjoint ^ 

Therefore, if all fields are in the adjoint of SU(N C ), one gets: 



^adjoi 



oint 



7l s 2 

— — 4n v + — lit 
6 3 i 



N c , (4.4) 



whereas if they are in the fundamental representation the corresponding value of b is: 

77 77 -f 

fundamental = ^ _ g + >JJ_ ^ 

12 v 3 V 7 

As discussed above, in the vector multiplet of (4,4) supersymmetry in two dimensions all 
fields are in the adjoint and n s = 4, rif = 2 and n v = 1. Therefore, one has: 

^vector multiplet _ _2j\T (4 5) 

On the other hand, for a hypermultiplet in the fundamental representation n s = 4 and = 2 
and the value of b in this case is given by: 

^hypermultiplet _ ^ (4 7) 

The total value of b in a Af = (4, 4) gauge theory with Nj matter hypermultiplets is thus 
—2N C + Nf and the running coupling constant is: 



1 1 



9ym 



9ym ( N Nl 

2yl c 2 



Eq. (14.81) shows that the (4,4) theory has negative beta function when Nf < 2N C , while for 
Nf > 2N C the beta function changes its sign and becomes positive. In the borderline case 
Nf = 2N C the one-loop beta function vanishes and the coupling does not run any more in 
perturbation theory. In the next subsection we will be able to reproduce the running (14.81) 
from a probe analysis of our gravity solutions, both from the unflavored and backreacted 
geometries, and, therefore, to characterize geometrically how the different regimes of iV c and 
Nf are encoded in our solutions. 

Some other interesting non-perturbative aspects of 2d M = (4, 4) theories, mostly related 
to conformal IR points, were studied in [32} [33], [31] . 
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4.1 Probe analysis 



In order to extract the information on the gauge theory contained in the gravity dual we 
have found, let us study the dynamics of a D3-brane probe moving under the influence of 
the metric and RR form of the background. The action of such a probe will be given by the 
standard sum of DBI and WZ terms: 



S 



-T, 



det(G 4 + 2na'F) + T 3 / G 



(4.9) 



where £ a (a = 0, ■ • • , 3) is a set of worldvolume coordinates, F is the field strength of 
the worldvolume gauge field and the hat over G4 and C4 denotes the pullback over the 
worldvolume of the D3-brane. The probe we are interested in is a color brane probe extended 
along the directions £ a = (x° , x 1 , 6 , <p) and located at a fixed value of the a coordinate at a 
constant value of all the remaining coordinates. We shall first consider the configuration in 
which the worldvolume gauge field vanishes. For such a configuration, the induced metric 
on the brane worldvolume is given by: 



H~i dx\ 1 + 



zH- 



{d9f + sin 2 # (l + o- 



m 2 H 



cot 2 #)(d0) 2 . (4.10) 



The determinant of this induced metric is: 



det Ga 



z sin 9 



„ m 2 H „ 
1 + a 2 —7— cot 2 9 . 



m?H v z 
Moreover, the pullback of the RR four-form potential (I2.12p is: 

zsin.8 



Ca 



m 2 H 



dx° A dx 1 A dB A 



(4.11) 



(4.12) 



If we now add the DBI and WZ contributions as in (j4.9p . we get: 



S, 



pot 



-T 3 



d x 



z sin 9 
m 2 H 



„ m 2 H „ „ 
1 + a 2 cot 2 - 1 



(4.13) 



Notice that the right-hand side of (14.131) is minus the static potential between the stack 
of N c color branes and the additional probe. This potential is in general non-vanishing, 
which means that there is a non-zero force between these branes. However, the potential 
does vanish for o = 0, which should be interpreted as the supersymmetric locus of the 
branes inside the Calabi-Yau space. Notice that the branes can be at any point in the M 4 
directions. Let us parameterize these flat directions in terms of four cartesian coordinates 
Z* (i = I,-- - ,4) as follows: 

(dZ 1 ) 2 . (4.14) 



dp 2 + p 2 dn 3 3 
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Next, let us assume that we are at the no-force point a = and let us allow the Z l scalars 
to depend only on the Minkowski coordinates = (a; ,^ 1 ). Moreover, let us assume that 
we switch on the worldvolume gauge field F in such a way that its only non-vanishing 
components are those along the unwrapped directions x M . The lagrangian density of the 
DBI term of the probe brane action is: 



DBI 



m?H 



sin 9 



(2irry') 2 

1 + [ ——^HF^F^ + H{d^f 



(4.15) 



where all functions are evaluated at a = 0. Let us also re-express the transverse coordinates 
Z l in terms of the scalar fields if 1 of the gauge theory living on the color brane as: 



2vraV . 



(4.16) 



In order to get the gravity expression of the Yang-Mills coupling we have to expand the 
lagrangian density (I4.15P (integrated over the angular variables (9, <fi)) up to quadratic terms 
in the fields and remember that passing from the abelian to the non-abelian theory amounts 
to substituting: 

FiafF »» _> tr[F, u F^} 
and similarly for the scalars </?\ By doing so, we get: 



2 ^ 



(4.17) 



C 



DBI 



(27T) 3 (« , ) 2 ^3 



quadratic 

The canonical kinetic term for the gauge field is: 

1 



trlF^F 



i 



(4.18) 



(4.19) 



where we have introduced a dependence of the Yang-Mills coupling on a renormalization 
energy scale /x. This energy scale should be related to the holographic coordinate p (see 
below) . 

By comparing with the expression (14.181) obtained by expanding the DBI action and taking 



into account that (27r) 3 (a') 2 T 3 
Mills coupling: 



l/g s , we arrive at the following gravity value of the Yang- 
1 z(p,a = 0) 



9ym 



m 2 g s 



(4.20) 



Taking into account the relation (I4.16P between the coordinates Z l and the scalar fields if 1 , 
which do not acquire anomalous dimensions, it is natural to assume that the energy scale p 
is related to the holographic coordinate p as follows: 



p = 2n a' p . 



(4.21) 
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Let us now substitute in (I4.20p the value of z(p, a = 0) as given by (13.201) for the general 
flavored solution. For concreteness we will assume that p > pq or, equivalently, that p > tuq. 
To obtain the corresponding result for p < itlq one should simply take Nf = in the 
equations that follow. After writing z(p, a = 0) in terms of p, we get the following YM 
coupling: 



gy M {p) m 2 g s 4nm 2 Q 2np 2 V c 2 
Let us next define the bare UV Yang-Mills coupling as: 

1 _ N f 

9ym m2 9s 4nm\ 

Then, we get the following running coupling: 



N c - -± . (4.22) 



(4.23) 



9ym{^) 9ym V 271 V 



1 - 



9ym < at N 



(N e --J-)). (4.24) 



By comparing (14.241) with the result obtained in the previous subsection for the one loop 
running coupling constant in perturbation theory (eq. (14.81) ) we see that they coincide 
exactly. Moreover, it follows from (14.181) that the action for the scalar fields takes the form: 

' tr[«9^<9V] • (4-25) 



9ym(v) 



This implies that the metric of the moduli space in the Coulomb branch of the theory is of 
the form: 



ds 



2 



(dp 2 + p 2 dVt 2 ^j , (4.26) 



M 9 YM ^) 

which is the result also obtained in the field theory [32]. Thus, we have succeeded in recover- 
ing from our gravity solutions the relevant perturbative information of the Coulomb branch 
of the gauge theory living on the branes. 

Let us now consider the possibility of placing the brane probe at a ^ 0. We have already 
shown that if a is constant and non-vanishing the energy of the vacuum configuration is 
non-zero, signalling the breaking of supersymmetry of this configuration. One can confirm 
this conclusion by studying the kappa-symmetric embeddings of probes (see appendix |C|) . 
Actually, one can argue that moving away from the point a = is dual to switching on a 
Fayet-Ilioupoulos parameter on the gauge theory side. The rigid displacement of the color 
brane in the direction of the a coordinate gives rise to a configuration of the probe that breaks 
the supersymmetry of the a = point. In order to recover this supersymmetry we must allow 
some non-trivial profile of the o coordinate which, in some cases, can be interpreted as the 
recombination of the color and flavor branes. We will explore this possibility in subsection 
14.21 in our description of the string dual to the Higgs branch of the theory. 
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4.2 The Higgs branch 



In this subsection we will present a holographic description of the Higgs branch of the d = 2, 
J\f = (4, 4) gauge theory. We begin by recalling that the Higgs branch is a phase of the 
theory in which the quark fields Q and Q acquire a non-vanishing expectation value. It is 
also interesting to point out that the gauge group that corresponds to our brane setup is, 
actually, G = U(N C ) « SU(N C ) x U(l). The presence of the U(l) factor in G allows us to 
turn on a Fayet-Ilioupoulos (FI) coupling on the lagrangian which, as we will argue below, 
will induce a non-zero VEV of the quark bilinear QQ and, thus, will force the system to 
enter into the Higgs branch^. To see how this can happen, let us recall [35] that the FI term 
can be recast as an extra term in the superpotential linear in the adjoint chiral superfield. 
Let $ denote the U(l) component of such a field. Notice that $ is naturally coupled to the 
matter fields Q and Q by means of the interaction QQ<&. Then, the superpotential W can 
be written as: 

iy = QQ$-r$ + --- , (4.27) 

where the dots refer to terms which do not contain $, and r is the FI coupling. By extremizing 
W with respect to $, one gets: 

QQ = r , (4.28) 

which shows that a non-vanishing value of the FI coupling implies a non-zero value of the 
bilinear Q Q, as claimed. 

To figure out how this mechanism can be implemented in our approach, it is interesting 
to recall the analogous Hanany-Witten brane setup in the type IIA theory (see [36] for a 
review). In this setup the (4,4) two-dimensional theories are engineered by considering a 
stack of N c D2-branes, with a compact dimension in their worldvolume. Moreover, the D2- 
branes are suspended between two parallel NS5-branes. The flavor branes in this setup are 
Nf D4-branes located between the two NS5-branes. If the latter are at different position in 
their transverse space, a FI term is induced in the gauge theory living in the D2-branes, with 
the FI coupling being the relative displacement between the two NS5-branes. When the two 
NS5-branes are misaligned in this way, the only possibility to preserve SUSY is by making 
the D2-brane end on the D4-branes, in such a way that both types of brane are recombined. 
Thus, we learn from this analysis that the Higgs mechanism is described as the reconnection 
of color and flavor branes which, in our d = 2, JV = (4, 4) system, is induced by a non-zero 
FI term. 

Coming back to our branes plus gravity setup, it is clear that, in order to find a holographic 
description of the Higgs branch, we should look for an embedding of a D3-brane which could 
be interpreted as representing a recombination of flavor and color branes. This is precisely 

7 Though in two-dimensional field theories an operator cannot have a VEV, we should think about VEVs in 
this section in a Born-Oppenheimer semiclassical approximation sense. The solutions describe a semiclassical 
state that is almost not changing in time. 
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the purpose of appendix £3 where we find a general class of D3-brane embeddings that 
preserve all the supersymmetries of the background. In these embeddings the D3-brane is 
extended along a holomorphic curve inside the CYi manifold. We will argue that one of 
these embeddings has the right properties to be considered as the string dual of the Higgs 
branch. Moreover, in subsection 15.21 we will analyze the meson spectrum in the Higgs branch 
in the probe approximation. Similar analysis in other brane setups can be found in ref. [37J. 

We want to find D3-brane embeddings that preserve all the supersymmetries of the back- 
ground and that interpolate between color and flavor branes. Recall that both types of 
branes are extended along two different two-cycles of the CY2 cone. Indeed, the color branes 
are extended along (9, 4>) (at a = 0), whereas the flavor branes of the Coulomb branch extend 
along [if), a), at fixed values of 9 and 0. To generalize these configurations it is natural to 
consider a flavor brane in which 9 and <fi are no longer constant. Let us thus consider an 
embedding for the flavor D3-brane in which the worldvolume coordinates are the same as 
in section [3l namely (x°, x 1 , if>, a), and the scalars 9 and (f> have a non-trivial dependence on 
the other coordinates of the CY2, namely: 

9 = 9(i),a), = 0(^,(7). (4.29) 

In order to determine the D3-brane embeddings of the form (14.291) which preserve the su- 
persymmetries of the background one has to study the kappa symmetry of the brane probe. 
This analysis is performed in detail in appendix The final result found in this appendix 
can be nicely recast in terms of the following two complex coordinates of the CY^. 

Ci = tan {^j , C2 = asmOe'^ . (4.30) 

It turns out that any holomorphic embedding of the type Ci = /(C2) solves the kappa 
symmetry equations and, thus, preserves the supersymmetry of the background. In principle 
any holomorphic function / will solve the kappa symmetry conditions. However, in order 
to make contact with the field theory analysis, it is rather natural to restrict ourselves to 
embeddings characterized by a polynomial equation of the type: 

£f = constant , (4.31) 

where the exponents p\ and P2 are constant integers. Notice that, in the Coulomb branch, 
color branes are extended along (1 (at ( 2 — 0), whereas the flavor branes span ( 2 at fixed 
Ci- For this reason, when pi = p 2 = 1 in (I4.3ip the corresponding embedding is symmetric 
with respect to the directions of the color and flavor branes that are recombined in the Higgs 
branch. Notice also the similarity between the embedding equation in this case, namely 
Ci C2 = constant, and the F-term equation (14.281) of the field theory. Indeed, the directions 
corresponding to the four scalars in the matter hypermultiplet can be identified, in our 
gravity dual, with those of the CY 2 , and then one can regard (14.311) for p\ = p 2 = 1 as the 
geometric realization of (14.281) in our holographic setup. 
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In the next subsection we will proceed to explore another interesting observable in a 
quantum field theory: the entanglement entropy. 



4.3 Entanglement entropy 

In quantum field theory the entanglement entropy between two complementary spatial re- 
gions A and B is defined as the entropy seen by an observer in A who does not have access 
to the degrees of freedom of B. It can be calculated from the reduced density matrix ob- 
tained after taking the trace over the degrees of freedom of B. In ref. [23] the authors 
proposed a simple geometric method to compute the entanglement entropy in quantum field 
theories which admit a gravity dual (see also [21]). This method consists in finding the 
eight-dimensional surface S with minimal area such that its boundary coincides with the 
boundary of A. Then, the entanglement entropy between A and its complementary region 
B is given by the integral: 



S = — / d^e- 2 ^detGs , (4.32) 

where G\q is the ten-dimensional Newton constant, given by G\q = 87r 6 c/ 4 g;? and Gs is the 
induced metric on E. We will consider a constant time surface S, obtained by minimizing S 
over all surfaces that approach the boundary of A at the boundary of the ten-dimensional 
bulk manifold and that are extended along the remaining spatial directions. 

For a gravity dual of a two-dimensional field theory the natural choice for the region A is 
just a line segment in the Minkowski direction x = x 1 . We will denote by / the length of this 
segment and we will choose our coordinate system such that A is given by {—1/2 < x < 1/2}. 
Moreover, for simplicity we will consider a surface £ that does not penetrate much in the 
IR region of the metric, in such a way that the UV metric (12.311) can be used. Furthermore, 
we will parameterize this 8d surface by the coordinates: 

C = (x,M,d,V,cO > ( 4 - 33 ) 
and we will assume that the surface is described by the function: 

u = u{x) . (4.34) 

It can be easily shown that this is a consistent ansatz for the minimal surface in the UV 
region of large u. By computing the induced metric Gg, we get that S is given by: 



vr 4 fk 



S = — - — / dxu 

m 6 G l0 



r2, m 4 

u H u 



(4.35) 
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Since the function S does not depend on x, the Euler-Lagrange equation derived from S can 
be integrated once and the result can be written as: 



u 



I 2 



u 



+ ^U 4 



u 



o • 



(4.36) 



where uq is the minimum value of u. For consistency we have to assume that uq is also large 
and we are always in the region in which the UV form of the metric (12.311) is valid. From 
the expression (14. 36[) we can obtain function of u: 



u> = ±^JUL)-l. (4.37) 

From this result we can compute the length I as a function of the turning point uq of the 
holographic coordinate. One has: 

l = 2 r* = h£r /- . ,4,8) 

It is now convenient to change the integration variable from u to £ = u/uq\ 

Z = M* f° J_ d £ . (4.39) 

This integral can be performed exactly and one can obtain the relation between I and Uq. 

_ 2^f^ r(§) i 



m 2 



r(|) «o ' 



(4.40) 



We can use (14. 37ft to eliminate vl in the entropy functional S. The resulting integral is 
divergent if the upper limit is infinite. Let us regulate this divergence by integrating up to 
some value Uoo of u. In terms of the variable £, one gets: 

ff = „e (44i) 

m"G 10 J 1 

The integral appearing on the right-hand side of (I4.4ip takes the value: 



r 4 ^=-f^)Vf-i,i,',f^) 6 )-^^f. (4.42) 



v /|e3i 2V Mo / V 3'2'3'V Moo ; ; 2 r A 
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Sending — > oo we obtain a term that diverges and terms that are finite: 



/ d^ 1 -^) 2 -y^ T -^ + o((^y) . (4.43) 

VFT 2 ^ U ° J 2 r (l) \ K uJ J 1 ; 



Plugging this result into (14.41 j) . we get the following equation for the entropy: 

4 4 
F 9 7T 
7/ — 

m 6 G 10 00 m6G 10 r(i) 



S ~ I^TT" M - - -TtV TTTTY^o ■ (4-44) 



As pointed out above, the first term in the above expression is divergent. Let us consider 
the finite part of S, obtained by subtracting the divergent part: 

Moreover, let us use f!4.40f) to write the value of S* mite in terms of the lenght I. After some 
calculation we get: 

r(|) 

where V* is the volume of the sphere of radius sfzlj m along which our D3-branes are wrapped: 



S finite {l) = -2v^F 



" 3 N 2 

V, -f , (4.46) 



V. = ^ . (4.47) 

Notice that the entropy (14. 46 ft scales as N 2 /l 2 . Actually, the expression we have just found 
is the one that corresponds to a 3+1 dimensional gauge theory compactified on a two-sphere 
(compare with eq. (42) of [24J ) . Indeed, from the fact that our solution is obtained from a 
stack of wrapped D3-branes, this UV behavior was to be expected. It would be interesting 
to explore the IR behavior of the entropy for our model. In particular one should investigate 
if, as happens in [21], there is a critical value of I such that the behavior (14.461) is modified. 
In our case finding the minimal surface in the IR region would require us to generalize the 
ansatz (14.341) for the embedding (the coordinate u would have to depend on the coordinate 
a defined in ( 12.271) ). We will not attempt to perform this generalization here. 



5 Meson spectrum 

In this section we discuss general aspects of the computation of the spectrum of mesons. 
By now there is a well known technique based on considering a (stable) probe brane in 
a given background and studying the normalizable fluctuations of the fields on the brane. 
The problem is of Sturm-Liouville type and typically admits a discrete set of normalizable 
fluctuations. 
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The string computation described above corresponds to the so called 'quenched approx- 
imation' in lattice field theory. Indeed, what happens in this case is that the background 
encodes the non-perturbative physics of the field theory with SU(N C ) gauge group, and the 
dynamics of the flavors (encoded on the 'probe' flavor branes) does not modify the color 
dynamics. This was discussed in many papers and it is known to correspond, in the dual 
field theory, to ignoring diagrams where fundamentals propagate inside the loops. The jus- 
tification for this approximation is that this type of diagrams are weighted by a factor of 
Nf/N c (see figure 3 and section 7 of [21J for a detailed discussion on this subject). It is 
observed in lattice simulations that the results of the quenched and unquenched calcula- 
tions of the masses of mesons do not differ qualitatively. Nevertheless, in some cases the 
numerical differences are relevant to match experimental data. This is one of the reasons 
why lattice people insist on more expensive unquenched computations. Besides, it is only 
for the meson spectrum that the quenched approximation works qualitatively well, while for 
thermodynamical, finite densite or other properties it gives the wrong physics. 

In this paper and in other papers along these lines, the backreaction of flavors has been 
considered. Hence it is natural to ask how to compute the mass of the mesons in this more 
realistic scenario. 

In principle, we can use what we know about computing the spectrum of glueballs. Fol- 
lowing the ideas introduced in [38] and elaborated on in [39] , we need to fluctuate the 
fields of type IIB supergravity and find sub-sets of fluctuations that constitute a consistent 
truncation. Remind that in general the glueball states may be linear combinations of the 
fluctuations of metric, dilaton and RR fields that diagonalize the linear fluctuation equations. 

In the case of backreacted flavors, we should find a consistent sub-set of fluctuations that 
includes closed and open string fields (the supergravity fields and the fields on the brane). 
More concretely, if hatted fields indicate the background solutions we discussed in the first 
sections of this paper and e a small parameter associated with the small fluctuation, we will 
have: 

Gfj, u — » G^v + edG^, F 5 — > F 5 + eSF 5 , 

C 2 — > e5C 2 , B 2 — > e5B 2 , 4> — ► , X — ► ^X , 

4 - eM« , if 1 -> , (5.1) 

where (A* , ip l ) denote the fields on the flavor branes with a flavor index i. So, the procedure 
should be to expand the action 

S = Sub + Sdbi+wz , 

up to order e 2 , find a subset of fluctuations that is closed and diagonalize them. It may 
happen that there is a mixing of glueballs with fields on the branes. Clearly the glueballs 
cannot have a flavor index, so the diagonal glueballs should be the same with or without 
flavors but the mesons can have a component of the 'closed string fields". 
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The procedure above is quite lengthy, but we believe that it should be the one to ap- 
ply to have a correct equation for the mesons and glueballs masses. We can nevertheless 
approximate this complicated computation by neglecting the mixing between glueballs and 
mesons. This approximation is similar to the one made when studying multi-electron atoms, 
and solving the Schrodinger equation of the external electron in the presence of the field 
created by the nucleus and the other electrons. Here, the "electrons" and "nucleus" are the 
flavor and color branes respectively §. This is equivalent to considering a probe brane in the 
background that encodes the dynamics of the field theory with gauge group SU(N C ) and 
flavor group SU(Nf). Doing this probe brane computation will indeed take into account 
the diagrams that scale like Nf/N c , because they are encoded in the background where the 
probe fluctuates. 

There is another interesting point: in the particular backgrounds we studied in the first 
part of this paper, the Nf quarks we added are massive, their mass being related to the 
radial position pq where we localize the flavor branes. This implies that if we place this 
probe brane at a fixed position p p above the value of pq the influence of the flavors will be 
explicit in the differential equations satisfied by the fluctuations, while if we place it below pQ 
there will be no influence whatsoever. This is indeed what happens in this and many other 
theories; one can choose a regularization where the matching between the high energy and 
the low energy theory is made exactly at the point pq. This seems to be the regularization 
preferred by our set-up. 



5.1 Quenched mesons on the Coulomb branch 

As argued above, when Nf « N c it is justified to neglect the backreaction and consider 
the flavor branes as probes in the background created by the stack of color branes. This 
defines the so-called quenched approximation, which has been used frequently to determine 
the spectra of mesons of the theory by analyzing the normalizable fluctuations of the flavor 
brane probe. In this subsection we shall illustrate how this can be done for our unflavored 
background. To begin with, notice that the metric Q a b induced on the worldvolume of a 
D3-brane embedded in the ten-dimensional metric ( 12. II) as in eqs. ( 13. 2p and ( 13. 3ft is: 

GatdCde = [H(p Q ,a)]- l > dxl x + ^ H ^ a) } 2 Ida 2 + a 2 (#) 2 l . (5.2) 

Let us now study the form of this metric when the worldvolume holographic coordinate 
cr — > oo. For this purpose we can use the asymptotic expressions of z and H obtained in 
section [2] (eqs. (I2.25P and (I2.26P ). One gets that, at leading order, the metric (15.21) reduces 
to: 

mV dx\ ! + — - (da) 2 + — (#) 2 , (5.3) 



3 We are grateful to Angel Paredes for suggesting this analogy. 
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where we have redefined the Minkowski coordinates ' . Notice that the 

metric ( 15.31) is of the form Ad S3 x S , with the S 1 factor being given by the ip coordinate. 
The radii of the AdS^ and S 1 factors in (15.31) are the same and equal to l/m. This UV 
behavior will be useful to establish the operator/fluctuation dictionary in our case. Notice 
also that, as H(pQ,a) and z(pQ,a) go to a finite value for a = 0, the circle parameterized 
by tp collapses in the IR at a = 0, as usually happens in the induced metrics that describe 
holographic flavor [10J. 

Let us now study the fluctuations of the flavor brane probe around the Coulomb branch 
configuration ( 13. 31) . To simplify matters we will only consider the fluctuations of the coordi- 
nate p around the point p = pq, namely: 

p = p Q + p{x»^,a) . (5.4) 

By expanding the DBI+WZ lagrangian density of the flavor brane up to second order in the 
fluctuation p , one gets: 



c = -Ti H l 2^^g g ab d a pd bP = 



(d x ,p) 2 + a(d a p) 2 + -(<V) 2 . (5.5) 



z a 



The equation derived from this lagrangian is given by: 

d a [ad a p] + ^-d 2 xl ,p + -dfyp = . (5.6) 

In order to find the solutions of this equation, let us separate variables in this equation as: 

p = x (a) e ikx e W , (5.7) 

where I is an integer, which for simplicity we will take to be non-negative, and A; is a 
momentum vector along the Minkowski directions. The equation resulting for x( a ) i s: 

d a [ad aX ] + \ — M 2 - -] X = 0, (5.8) 
L J L z u \ 

where M 2 = —k 2 . Our purpose is to determine the values of M for which there exist 
regular solutions of (15.81) . Notice that a is a non-compact variable that can be arbitrary 
large and, when pq is large enough, its minimal value is o = 0. Moreover, from the behavior 
of the functions z(pq, a) and H(pq, a) for a = 0, 00 one can verify that, in both limits, the 
fluctuation \ behaves as a linear combination of a 1 and o~ l for I 7^ 0. When I = the two 
independent solutions behave as \ ~ (constant, log a) as o —>■ 0, 00. Therefore, it is clear 
that, by fine-tuning M, there are solutions of ( 15.81) that are regular when a = 0, 00. The 
discrete set of values for which this occurs is identified, in the context of the extension of the 
gauge/gravity correspondence put forward in |10j . with the mass spectrum of the mesons 
in the dual field theory. Actually, since the induced UV metric is AdS% x S 1 , one can use 
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the full machinery of the AdS/CFT correspondence to determine the conformal dimension 
of the operator dual to the normalizable fluctuations. From the behavior at large a of the 
solutions of (15. 8p one easily concludes that the p fluctuations with angular quantum number 
/ are dual to operators with dimension A = I + 1. 

The values of the meson masses can be obtained numerically by means of the shooting 
technique. In tabled] we present the values found by this method for I = 0. 



Meson spectra for I = 


n 


WKB 


Numerical 


1 


4.23 


4.22 


2 


7.32 


7.31 


3 


10.35 


10.35 


4 


13.37 


13.36 


5 


16.37 


16.36 



Table 1: Numerical and WKB values of M for 1 = 0. These values correspond to pq = 1.5 
and = 5. 

It is also interesting to point out that the fluctuation equation (15.81) can be converted into 
a Schrodinger equation by means of a suitable change of the holographic variable. Indeed, 
let us define the variable y as follows: 

e y = a . (5.9) 

Notice that, in this change of variables a — > oo is mapped into y — > oo, while a = 
corresponds to y = — oo. In terms of y, the equation (15.8]) of the fluctuations can be written 
as the zero-energy Schrodinger equation: 

0-^(y)x = o, (5.io) 

where the potential V(y) is given by: 

V(y) = l 2 - M 2 - — — . (5.11) 

z 

It is easy to verify by using the behavior of the functions H and z at a = 0, oo that V(y) — > I 2 
when y —>■ ±oo. In figure [3] we have plotted the form of this potential for I = 1. Notice 
that V has a unique minimum and that, for 1^0, the classically allowed region V < of 
(I5.10P is of finite size. Then, it is clear that the wave equation for V has a discrete spectrum 
of energies and, for certain values of M one of these levels has zero energy as in (15.101) . A 
fast, and rather accurate, estimate of the mass spectrum can be obtained by applying the 
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Figure 3: The Schrodinger potential for / = 1 and M 2 = 18.11, which is the mass of the 
ground state (n = 0) for this value of I. We have set = 5 and pq = 1.5. 



semiclassical WKB method. This method is explained in detail in [lOj HI] (and applied to 
compute meson masses in [42, 43J). Here we will limit ourselves to give the final result for 
the WKB spectrum, which is: 

M 2 WKB = f 2 (n + l)(n + 21) , (5.12) 

[C(pq)\ 



where n is a non-negative integer and C{Pq) is the following integral: 



H(p Q ,a 



C( PQ ) = / da J ™ / . (5.13) 



a 



In (I5.12p n > when 1^0, whereas n > 1 for 1 = 0. Let us now define the mass gap M* 
of the mesonic spectrum as the mass of the lightest meson. We can get an estimate of this 
mass gap from the WKB formula for the masses. In this case the lightest meson corresponds 
to taking I = and n = 1 in (15. 12ft . One has: 

MY KB = ^ ■ (5.14) 

We can further estimate the value of the integral C(Pq) i n (15.131) by using our asymptotic 
expressions (12.251) and (j2.26|) : 



rn 

This integral can be explicitly evaluated, namely: 



30 J 

C(P Q ) « ^ / 2 Z 5 • (5.15) 



CM « ^ ■ (5.16) 
30 



Notice that it is independent of z*. Using the value of C{Pq) obtained above, we find: 

2 v / 2m Q 



M 



WKB 



2V2- 2 



m p Q 



(5.17) 



aA g s N c 

Thus, we predict that the mass gap is linear in the quark mass iriq and is independent of the 




Figure 4: Numerical mass gap as a function of pq. The straight line corresponds to the 
approximate WKB prediction. We have set = 5. 



IR parameter z*. Notice that, in principle, the use of the asymptotic expressions (I2.25p and 
( I2.26P is only justified when pq is large. However, it turns out that (15.171) is a rather good 
approximation to the actual numerical results (see figure H]). Also, the numerical results show 
that the mass gap is almost independent of z*, as predicted by our approximate analysis. 
Furthermore, by using (15.161) in (15.121) we also obtain an approximate analytical expression 
of the WKB meson mass spectra, namely: 



M 



WKB 



(n,l) « 2m 2 p Q y/(n + l)(n + 2l) . 



(5.1* 



This formula works rather well for I = 0, 1. 

Up to now we have only considered one particular fluctuation mode of the flavor brane. 
The remaining modes can be treated similarly and one can find decoupled sets of fluctuations, 
for which one can obtain the corresponding mass levels. Let us illustrate this fact with a 
particular example. Let us consider the fluctuation of the coordinates 9 and that determine 
the position of the flavor brane in the Calabi-Yau manifold, namely: 



e = e Q + 6{x»,iI),<t) 



= 0o + (f)(x fl ,ip,a) 



(5.19) 



where #0 an d 4>o are the constant unperturbed values. By analyzing the behavior of the 
DBI+WZ action for these perturbations one discovers that these fluctuations are coupled. 
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However, following the same steps as in [15] . one can decouple them. Indeed, let us define 
X+{ a ) and X-{ a ) i n terms of 9 and by means of the following relation: 

= \ (x+(<x) + X-W) sin O sin(^) , 

0=1 - X ») e ik * cos(ty) , (5.20) 

where I is a non-negative integer. By explicit calculation one can verify that x+ an d X- 
satisfy decoupled equations that, actually, can be mapped to ( 15. 8ft by a simple redefinition. 
Moreover, from the AdS/CFT dictionary one can associate to x+ a field of conformal di- 
mension A + =1 + 3, while X- is dual to a field with dimension A_ = 1 — 1. It is interesting 
to notice the similarity between these results and the ones found in [131 HI] for the D3-D3 
intersection in flat space, which is dual to a 2d defect field theory. 



5.2 Meson spectra on the Higgs branch 



Let us now consider the fluctuations of the brane probe around a non-trivial holomorphic 
embedding of the type ( 14.31 j) . As in the analysis of the Coulomb branch in subsection 15.11 
we will concentrate on studying the fluctuations of the coordinate p around p = pq. If p 
denotes this fluctuation (see eq. (15.41) ). the corresponding lagrangian up to second order in 
p takes the form: 

T 3 



C = 

2 



Ax (d x »p) 2 + A a (d a p) 2 + Aip {d^pf 
where the coefficients A x , A a and A^ are given by: 

mV 2 (1 + p cos 6) 2 H + p 2 z 2 sin 2 8 



(5.21) 



Aip 



m 2 a (1 + pcosl 
a (1 + pcos9) , 



(5.22) 



with p = pijpx and 9 being the function of a given in (1C.18I) and H and z are evaluated 
at p = pq. As it should, when p = the lagrangian (I5.2ip reduces to (|5.5p . Moreover, the 
equation of motion derived from (I5.2ip is: 



A a d a p + A x d 2 ,p + — d%p = 



(5.23) 



From now on we shall restrict ourselves to the analysis of the fluctuations around the em- 
bedding with p = 1 which, as argued above, we think realizes geometrically the brane 
recombination of the Higgs branch. In this p = 1 case, the function 9(a) is given by: 



sin 



a* 
a 



(5.24) 
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where cr* is a constant that represents the minimal value of the coordinate cr (that occurs for 
9 — 7i"). Notice that, when a* 7^ 0, the probe brane never reaches the point a = 0. In this 
case the brane starts at the north pole 9 = when a —>■ 00 and as a — > cr* it approaches the 
south pole of the (#, 0) two-sphere. From the point of view of the color brane that is being 
recombined the fact that o = is not reached means that a FI term has been switched on. 
According to our field theory discussion this triggers the Higgs mechanism. 

By using (15.241) in ( 15.221) the coefficients A x and A a for this p = 1 case can be easily 
written in terms of the coordinate cr, with the result: 



A x 



2(o-o*)H 2o*z 



m 2 cr 2 



A a 



2 (a - a* 



and the differential equation for the fluctuations becomes: 



(cr - cr*) d a p 



+ 



(a-a^H a^z 



m 2 a 2 



d l»P + 77 7 d %P = 

4(cr - a*) v 



(5.25) 



(5.26) 



After separating variables as in (15. 7p . this equation becomes: 



(a - o-*) d a x 



cr 



cr* 



H o*z 
— + 



m 2 a 2 



M - 



4(cr 



X = . 



(5.27) 



Eq. (15.271) should be compared with the analogous equation (15. 8p for the Coulomb branch. 
Curiously, when a* vanishes eq. (I5.27P becomes (15. 8p with / replaced by 1/2. This makes 
sense because 9 is constant (and equal to zero) when a* = is taken in (15.241) . However, the 
embedding for cr* = is not the same as the one in the Coulomb branch because the angle 
<p is not constant (see (IC.16P for p = 1). Somehow, a* = is a singular point in the family 
of p = 1 embeddings we are considering. Clearly, for this singular point, the mass spectrum 
of the Higgs branch is related to the one in the Coulomb branch by means of the following 
relation: 

M Hi " s (n,l) a , =0 = M Coulomb (n, ~) , I e 21 . (5.28) 



Notice that the spectral flow relation (I5.28P only makes sense if / is even. The numerical 
calculation of the mass spectra by using the shooting technique shows that the spectrum for 
a* very small differs significantly from the one corresponding to cr* = and, therefore, the 
mass levels are discontinuous at a* = 0. Actually, when a* is increased from values very 
close to zero, the masses vary rapidly until they reach a plateau where they become nearly 
independent of a*. We have illustrated this fact in figure [5j In order to have a clearer idea 
on how the masses behave when a* is not very small, let us apply the WKB method also in 
this Higgs branch case. With this purpose we first map (15.271) into a Schrodinger equation. 
First, we define the new variable y as follows: 



cr — cr* 



— 00 < y < +00 . 



(5.29) 
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I t t tl 



♦ 1 = 

. /= 1 



log I0 (o»m) 



Figure 5: The mass of the ground state in the Higgs branch as a function of a* for I = and 
1 = 1. The straight line corresponds to the prediction of eq. ( 15.351) . We have set pq = 1.5 
and = 5. 



Let us also define y* as: 



(7* 



(5.30) 



Then, one can demonstrate that (15.271) can be written as the Schrodinger equation (15.10!) . 
with the potential V(y) being given by: 



V(y) 



4 



z 



e 2y H + 



?y-y* 



(l + ev-y*) 



(5.31) 



As a check, notice that when y* — > — oo (which corresponds to taking a* — > 0) the potential 
(15.31 j) becomes (15. lip with 1^1/2. Interestingly, the mass spectrum associated with the 
potential (I5.3ip has exactly the same form as in (15.121) with a different value of the prefactor 
C, which is now given by the integral: 



C(PQ,cr*) 



da 



(T,, 



m 2 a 2 (a — (X* 



(5.32) 



Let us now estimate the mass gap in the Higgs branch. We first notice that, if <r* is not very 
small, the second term inside the square root in eq. ( 15.321) dominates the integral. Thus: 



m 



da 



^z(p Ql a) 
a\f< 



a — a* 



(5.33) 



Taking the further approximation in which z is taken to be constant, which will be valid for 
large p Q , we get: 



m 



da — -j== 

a^Ja-a* 



TC./Z 



m 



(5.34) 
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which is independent of pq. Using this value we get that the mass gap in the Higgs branch 
(for large pq) is approximately given by: 




(5.35) 



5 10 15 20 



Figure 6: The mass gap in the Higgs branch as a function of z*. The continuous curve 
corresponds to the prediction of eq. (I5.35p . We have set pq — 1.5 and a* = 10. 

Notice that in this case the gap is independent of the quark mass and only depends on 
the infrared scale Moreover, in the decompactification limit — > oo the mass gap 

vanishes, in agreement with the result found for the co dimension- two defects in the D3-D3 
intersection in flat space (where the discrete spectrum of mesons is lost) [37]. In figure [6] we 
compare the numerical values for the mass gap with those obtained by means of eq. (I5.35p . 
We see that the agreement between our approximate estimate and the numerical values is 
excellent. 



6 Additional supergravity solutions 

Let us now try to find new simpler solutions of the BPS equations (12.91) . We shall try an 
ansatz in which the three functions g, H and z depend on the two variables p and o as 
follows: 

g = g(a) , H = H(p) , z = Zl (p) z 2 {a) . (6.1) 

Notice that, with this ansatz, the third of the BPS equations (12.91) is automatically satisfied 
(both sides of the equation vanish). Moreover, by plugging the ansatz on the first equation 
in (12.91) and by grouping in the same side the terms that depend on the same variable, we 

^ = -2c,, (6.2) 
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where we have called the separation constant — 2c\. Eq. (16.21) contains two different equa- 
tions. The one for z\{p) is the following differential equation: 

z[(p) = , (6.3) 
P 6 

which can be integrated as: 

zi{p) = c 2 + % , (6.4) 
P 2 

with c 2 being a new constant. We also get from (16. 2p the form of g(a) in terms of z 2 (p), 
namely: 

g(a) = -^z 2 (a). (6.5) 

From the second of the BPS equations ( 12. 91) . we obtain: 

m 2 H(p) z 2 (a)z 2 (a) 
z((p) a 

where C3 is another constant. The resulting equation for z 2 (a): 

z 2 (o)z 2 (o) = c 3 cr , (6.7) 

can be integrated as: 

z 2 (a) = a/q + c 3 a 2 . (6.8) 

Moreover, from (16.61) we also learn that the warp factor H(p) is related to z\(p) as follows: 

H(p) = — 2 zl{p). (6.9) 

Using the value of zi(p) and z 2 (p) given in eqs. (16.41) and (16.81) . we find that the solution we 
have is: 

2ci 



9(0) = o V c 4 + c 3 a 2 



H{p) 



C 3 r ^ t2 



Ci 

C 2 + 



z(p, a) = (c 2 + j 2 ) + c 3 a 2 . (6.10) 

The ten-dimensional metric corresponding to this solution can be obtained by plugging the 
values of H and z given in (16.101) into the general ansatz (I2.ip . One gets: 

ds\ Q = dsl + ds 2 , (6-11) 
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where the six- dimensional part of the metric ds\ is independent of a and given by: 

ds\ = dxl A + H s dp 2 + p 2 dtt 2 3 
and the four-dimensional metric, which does not depends of p, can be written as: 

C3 



(6.12) 



2 a/q + c 3 a 2 
ds A = 



«V C 3 



dO z + sin^ « 



+ 



m 



c 3 cr z 



At 2 + cr 2 ( dip+cos 6d(p ) ' 



(6.13) 

In order to rewrite the metric (16.131) in a more familiar form, let us perform a change of 
variables and define a new radial variable (, related to a as follows: 



c 2 



\/c 4 + c 3 a 2 . 



«V C 3 

Actually, if we define the new constant a as follows: 

_ 4 16c 4 

9 ' 

m z C3 



a 



the relation that gives a in terms of ( is the following: 



m r— T 

o- = T VC 4 - a 4 



(6.14) 



(6.15) 



(6.16) 



Clearly, £ > a, which corresponds to the range a > 0. After this change of variable, one can 
easily prove that the metric (16.131) becomes: 



ds\ 



d( 2 



(!) 



^ + 4 



d9 2 + sin 2 #cM 2 + 1 



c 



(# + cosfl#) 2 , (6-17) 



which is the metric of an Eguchi-Hanson space EH± with resolution parameter a. When 
a = the metric ( 16.171) becomes the one corresponding to the C 2 /Z2 orbifold. Moreover, if 
we define the new constants rj and R as: 



m 



then, the warp factor H becomes: 



H 



R 2 



C3C1 



m 



Using this result the six-dimensional part of the metric takes the form: 

doc 1 



dsl 



V 



+ 



dp 2 + p z dQ 2 3 



(6.18) 



(6.19) 



(6.20) 
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Notice that by rescaling x 0,1 and p in (16.201) we can have r\ = 0,1. The case r] = corresponds 
to a ten-dimensional metric of the form AdSs x S 3 x EH4, with R being the common radius 
of the AdSs and S 3 parts. When rj 7^ the six-dimensional metric flows in the UV region 
p — > 00 to the flat six- dimensional Minkowski space, whereas in the near-horizon region 
p — > it becomes AdS% x S 3 . Notice also that the functions g and z can be written in terms 
of the new variables as: 



9(0 = -f C 2 



*(p,C) 



( R \ 2 



4 



(6.21) 



Let us now write the expression of the RR five-form F 5 for this solution. In order to do that, 
let us choose the following vierbein one- form basis for the EH4 metric (16. 17ft : 



£ d9 
2 



- sin#< 
2 



a\ 4 



( dip + cos 9 d(f>) , 



d( 



(?)' 



(6.22) 



In terms of this basis, the Kahler form of the EH4 space takes the form: 

J = rj 1 A rj 2 + rj 3 A rf . 



(6.23) 



To get the expression of F 5 we just plug the value of g for our solution, displayed in ( 16.101) . in 
our ansatz of eq. (12.71) . After writing the result in terms of the variable ( defined in (I6.14p . 
one gets: 

F 5 = 2R 2 ( —- dp A dx° A dx 1 - u 3 ) Aj, (6.24) 
V p 6 H / 

where J is the two- form defined in (16.231) . Interestingly, by using the explicit value of H 
written in (I6.19p . one can verify that, in the AdS% x S 3 limit rj = 0, _F 5 takes the form: 



2 
R 



Vo\{AdS 3 ) - Vol(S 3 



A J , 



(6.25) 



which seems to be the same type of solution obtained in 



6.1 Non-relativistic backgrounds 

According to the papers in [16], one can construct a background dual to a non-relativistic 
theory in d + 1 dimensions by performing a combination of suitable T-dualities and shifts 
to a solution of supergravity in d + 3 dimensions. Here we shall apply this procedure to 
our unflavored background. The final result will have the isometries of the Galileo group in 
d = and, thus, it is somewhat dual to a matrix model. 
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As the first step in our procedure, let us introduce the light-cone variables x ± as: 



x ± = x° ± x 1 



(6.26) 



Next, we will perform a T-duality along the fiber direction d^p, followed by a shift in the 
light-cone coordinate x~ of the form: 



dx — > dx + •ydifj . 



(6.27) 



As the result of these two consecutive operations we obtain a background of type IIA super- 
gravity. The corresponding metric and dilaton are given by: 



dsj IA = H 5 



— dx + dx~ — r ydx + dip + — - ( d6 2 + sin 2 6 dip 



+ 



■ da A + dp z + p 2 dttj 



+ H 



.1 z 



di? , 



2* 



,2<I>o 



while the non-vanishing NSNS and RR potentials are: 
B 2 = cos 6 d(j) A dip , 



(6.28) 



C 3 = g 0J3 + — dx + A dx A da — 7 — dx + A da A dip , 
2z 2z 



a = 



2m 2 H 



dx + Adx A uj 2 A dip . 



(6.29) 



Let us next perform another T-duality along d^. We get a background of the type IIB theory 
with constant dilaton and the following metric: 



ds 2 IIB = H 5 



- dx + dx~ - ^ — (dx + ) 2 + —Add 2 + sin 2 6 d<p 2 ) 1 + 
4 z m z \ /J 

i (^da 2 + a 2 ^ + cos9d^y^j + dp 2 + p 2 dQ 2 3 j , (6.30) 



which has the following potentials: 



B 2 = --— dx + A (dip + cos 6d<p ) , 
2 z 



C 2 = -^-dx + Ada , 

2 z 



C 4 = g lo 3 A ( dip + cos #g?0 ) + — da A (dip + cos 9d<j)) 



2m 2 # 



UJ 2 



A dx + A dx~ 



(6.31) 
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Notice that, if we take the shift parameter 7 to be zero, the two-forms vanish and we recover 
the original ansatz we started with. From these potentials one can compute the corresponding 
field strengths as if 3 = dB 2 , F 3 = dC-i and F 5 = dC± + C 2 A H 3 . By computing directly the 
exterior derivatives and after using the BPS equations, one gets: 



H 3 = 1 - 
3 2 



m 2 a 2 g , 2a / m 2 a 2 H , 
dp — — I 1 — — da 



p 3 z 2 z V 2z 



T a 2 

A dx + A (dip + cos 9d(f) ) — dx + A u. 

' 2 z 



7 m 2 ag 



F 3 = 5—5- dx A da A dp , 

2 

F 5 = ^1 + A u; 3 A ( dtp + cos 6><i0 ) + g uo 2 A uj 3 . (6. 

As a check of the correctness of the final result one can verify that the forms satisfy their 
equations of motion, namely: 



dF R 



-H 3 A F 3 , d*F 3 = H 3 A F 5 , d(V 2 ^ *# 3 ) = F 5 A F 3 . (6.33) 



Let us interpret the effect of the transformations performed above, on our field theory, 
following what is explained in the papers [46J. 

As we mentioned above, if we take 7 = in ( 16. 27ft we recover the original background, but 
we impose that the coordinate x~ is periodic and hence it has discrete modding. Indeed, this 
would reproduce in string theory the "mass operator" of the Galileo invariant field theory 
that has discrete eigenvalues. In consequence, for 7 = we would be doing DLCQ of the 2d 
theory. Now, we can think of the nonzero 7 case as a particular case of a TsT transformation 
introduced in [J7]. In this particular example, it is generating a dipole field theory [48J. 
In a dipole theory we change the product between two generic operators $i(a;~,x + ) and 
$ 2 (x~,a; + ), with charges (gi,^), to the "star product": 

i( 8 - - d - ) 

$1(0:7, X+) * $ 2 (X2,xt) = e ^ 91 * 2 " ®l(xI,xt)®2(x2, X t) ■ ( 6 - 34 ) 

In our case (51,52) are R-symmetry charges (represented by traslations in the angle ip). As 
a consequence of this change of product, there will be physical procesess that are affected 
and some that are not. 

In particular, all products of operators that contain fields charged under the R-symmetry 
will change after the procedure described above. For example, masses of R-charged fields 
will change. All this will leave us with a non-relativistic field theory, that in the particular 
case we studied here, will also be non-supersymmetric (fermions, for example, will get a mass 
different from scalars and gauge fields will remain massless). 

As a word of caution, if we want to learn about non-perturbative aspects of this field theory, 
computing with the background in (16.321) . we should be careful to use the background for 
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values of the coordinates where the curvature is not very large. Note that since the x~ circle 
has zero lenght, strings that wrap the cycle are massless and the supergravity approximation 
may not be applicable. 

7 Fractional D1/D5 system 

In this section we will show how to describe the (4,4) two-dimensional gauge theory by means 
of a system of fractional Dl/D5-branes on the orbifold 1R 1,5 x C 2 /Z 2 . Fractional branes are 
a kind of brane that are specific to orbifold (and conifold) singular backgrounds. As can 
be checked by explicit calculations of string scattering amplitudes [19] , these Dp-branes can 
be viewed as ordinary D(p + 2)-branes wrapped on an integer basis of vanishing two cycles 
[50] . As in the analogous constructions of [51], this fact will be of fundamental importance 
to find the supergravity solution corresponding to the D1-D5 system we are investigating. If 
we parameterize C 2 by four cartesian coordinates x 6 ,x 7 ,x s ,x 9 , then the action of Z 2 will be 
simply given by: 

(x 6 ,x 7 ,x 8 ,x 9 ) — (-x 6 ,-x 7 ,-x 8 ,-x 9 ) . (7.1) 

Moreover, if the coordinates of R 1 ' 5 , the fractional Dl-branes are extended 

along a; , x 1 and are stuck at the orbifold fixed point x 6 = x 7 = x 8 = x 9 = 0. In this setup 
the flavor branes are D5-branes extended alone; X * X j els well as along the four directions of 
C 2 /Z 2 . Therefore, the full D1/D5 configuration is described by the array: 



c 2 /z 2 








1 


2 


3 


4 


5 


6 


7 


8 


9 


N c Dl 






















Nf D5 























To find a supergravity solution for this setup it is interesting to regard the C 2 /Z 2 orbifold as 
the limit of an ALE space (or, conversely, the ALE space can be thought of as the blow-up 
of the orbifold). This ALE space is endowed with an ant i- self- dual two-form u) 2 , which is 
the Poincare dual to the exceptional two-cycle E 2 that is obtained by the resolution of the 
orbifold singularity (S 2 shrinks in the orbifold limit). 

The peculiarity of the orbifold backgrounds is the appearance of six- dimensional twisted 
fields. They are p-form fields that can be seen as coming from higher dimensional (p + 2)- 
forms which have components along the non-trivial shrinking two-cycle S 2 of the orbifold 
at the orbifold fixed point. As such, they are the zero modes of the KK reduction of those 
(p+2)-form fields which have components along the anti-self-dual two-form cD 2 of the orbifold 
[52] . In our case the relevant twisted fields are the two- form RR potential A 2 and the NS-NS 
scalar b. They are related to the ordinary RR four-form potential C4 and to the NS-NS B 2 
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field as: 

C 4 = A 2 Alj 2 , B 2 = bCj 2 , (7.2) 
where the two form uj 2 is normalized as follows: 

Cj 2 A Cj 2 = — 1 . (7.3) 

C2/Z 2 

After performing this KK reduction, the action of the type IIB supergravity theory in the 
Einstein frame can be consistently truncated to be: 



S IIB = —^-\ / d 10 xV-det GR - - d$A*d$ - - I e*F 3 A* F 3 - 

^ K orb J ^ J 2 



- i / e~*d& A* 6 db-] [ T 3 A* 6 T 3 + ~ / A 2 A rf6 A F 3 ), (7.4) 

where we have not included the contribution of the different sources. In (17.41) the integrals 
in the second line are six-dimensional, = v2«4 = (2n)? g s (a') 2 , F 3 = dC 2 is the usual 
RR three-form field and T 3 is the following six-dimensional anti-self-dual three-form: 

T 3 = dA 2 - C 2 A db . (7.5) 

The anti-self-dualit}|^| of T 3 has to be imposed by hand. Actually, we will represent T 3 as: 

T 3 = T 3 - * 6 T 3 , (7.6) 

with T 3 being given by: 

T 3 = dA 2 - C 2 A db , (7.7) 

and d* 6 T 3 = 0. From the action (17.41) one can readily get the following equations of motion 
for the dilaton and forms: 

d*d<$> - ^e*F 3 A* F 3 + ]e~*db A* 6 db A Q 4 = , 

ci(e* *F 3 ) = db A T 3 A fi 4 , 
d(e-** & db) = T 3 AF 3 , 
d* 6 T 3 = -dT 3 = F 3 A db , (7.8) 

where S7 4 is the following four-form: 

Q 4 = 5{x 6 ) ■ ■ ■ 5{x 9 ) dx 6 A ■ ■ ■ A dx 9 . (7.9) 



J Our conventions are such that, in flat indices, e u " a = e u ■ " = 1 
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The Einstein equations derived from (I7.4p read: 



Rmn - ^G MN R = g (d M $d N <$> - ^G M Ndp$d p $ ) + 



+ 5 4 ^==£ |V* (d M bd N b - ^G MN d P bd p b^j + \tm MPQ T$ n \ , (7.10) 

where Gq denotes the restriction of the ten-dimensional metric to the six-dimensional sub- 
space R 1,5 , <5 4 = 5(x 6 ) ■ ■ ■ 5(x 9 ) and M, N = 0, • ■ ■ , 9. It is understood that in the terms 
where a 5<± appears in (17.101) (in the last line) these indices run just along the six directions 
transverse to the orbifold. 

The action (17.41) must be supplemented with the boundary actions corresponding to the 
different branes present in the problem, which give rise to source terms that must be added 
to the equations of motion (17. 8p and (17.101) of the bulk fields. In our particular problem the 
branes sourcing the bulk fields are N c fractional Dl-branes and Nf D5-branes that wrap the 
entire orbifold. Therefore, the boundary action is: 

S b = N c S m + N f S D5 . (7.11) 

In order to write the expression of Sdi an d Sdb, let us recall [52J that the NS-NS b field has 
a non- vanishing constant value for the C 2 /^2 orbifold. With our conventions, we can write 
b as: 

b = -(2ttV + 6) , (7.12) 

where b represents the fluctuation around the background value b = —2-K 2 a'. The action of 
a fractional Dl-brane in the Einstein frame is: 




Sdi = - — \ I d z xe--^-det(G af3 + e~ a 2na'F a(3 ) \ 1 + -^-j 




2n' 2 a' I 2ttV 



A 2 }. (7.13) 



Similarly, the action for the D5-brane is given by: 



+ y i(27v) { / d2xe ~" ^~ det ^ ~ b ~ J A *} ' ( 7 - 14 ) 
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where we have not included higher order terms in the twisted sector. In eqs. (17.131) and 
(I7.14p the tensions 7\ and T 5 are given by: 

^1 = o — / — ' = /o \5 / A3 — • (7-15) 

2na' g s (2n) b (a 1 ) 6 g s 



7.1 The unflavored solution 

First of all, let us now look for a solution corresponding to a stack of N c fractional Dl-branes 
on the space M 1,5 x C 2 /Z 2 , without D5-branes. With this purpose, let us make the following 
ansatd 10 ! for the metric and the dilaton: 

H- 3/ % aP dx a dxP + Hl'^dr 2 + j J>V + da 2 + jX)^') ' 

i=i j=i 

H{ /2 , (7.16) 

where r is the radial direction transverse to both the Dl-branes and the orbifold space, u l 
are the usual left-invariant one forms parametrizing the corresponding S 3 , a is the radial 
orbifold direction and u) 1 , with the range of angles conveniently chosen, are the left-invariant 
one-forms parametrizing the orbifolded S 3 . The left-invariant one- forms are normalized such 
that duj 1 = | eijkUji A u k and duj 1 = | e^cD- 7 A £j k . We will assume that there is a non-zero b 
field, that only depends on the coordinate r. Moreover, the warp factor Hi will depend on 
both radial coordinates. Then: 



H 1 = H 1 (r,a), b = b(r) . (7.17) 

In addition, our background will be endowed with a non-zero RR two-form C*2, as well as 
with a two-form potential A2 in the RR twisted sector: 

C 2 = (Hi 1 - 1) dx° A dx 1 , A 2 = a(r)dx° A dx l , (7.18) 

where a(r) is a function to be determined by the equations of motion. By using eqs. (17.61) 
and (17. 7p it is straightforward to find the field strength T 3 corresponding to the two-form 
potentials (17.181) . One finds: 

T 3 = [a' - b\H- 1 - 1)1 dr A dx° A dx 1 + [a' - b'{H^ - 1)1 u 1 A u 2 A to 3 , (7.19) 

8 

where the prime denotes the derivative with respect to r. Similarly the dot will denote the 
derivative with respect to a. 

10 a, (3 = 0, 1 label the worldvolume directions of the Dl-brane. 
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Let us now analyze the different equations of motion for this system. The equation of 
motion for F3 leads to the following equation for H\ : 

d r (r 3 a 3 H[) + ^(rVifx) = b'H ir 3 [a' - b'(H^ 1 - 1)] 5 4 - N^T^ , (7.20) 

where 5§ = S(x 2 ) • ■ ■ 5(x 9 ) and we have included the source term due to the presence of 
fractional Dl-branes. For the NS-NS twisted scalar b we get the equations: 

b' = -a' , (7.21) 



(b'r 



2vrV 



(7.22) 



where again sources have been included and £4 = 5(x 2 )5(x 3 )S(x 4 )5(x 5 ). Using this equation 
we can rewrite (17.201) as: 

d r (r 3 a 3 H[) + d CT (rV#0 + b' 2 r 3 5 A = -N c K 2 orh T x 5 8 , (7.23) 

One can verify that (17.211) . (17.221) and (17.231) imply the fulfillment of the equation of motion 
for the dilaton, as well as the Einstein equations (17. 10p (with the aditional source contribution 
coming from the boundary action). 

The solution of (17.221) is easily determined to be: 

Nc^T, 1 



b = -2ttV + 



2-n 2 a' 4n 2 r 2 



(7.24) 



where, to fix the constant term in (I7.24p . we have already taken into account that the 



background value of b is — 2ir 2 a f . Moreover, since K 2 rb Ti = (27r) 
(ELU as: 



la 



/\3 



-2-Ka' 



Ag s a'N c 



g s we can rewrite 



(7.25) 



To understand how our solution captures the gauge dynamics on the world- volume of the 
Dl-brane we will make a probe computation, similar to the one performed in subsection 14. II 
Accordingly, we shall consider a fractional Dl-brane probe in the background just described, 
in which a worldvolume gauge field F Q/ 3 has been switched on. The Dl-brane will be extended 
along sc ' 1 , at fixed values of the remaining coordinates. By expanding the action (I7.13P in 
powers of F a p, we find an expression of the form: 



S 



Dl 



-V 



^9ym 



d 2 xF af3 F a P + 



(7.26) 



where V can be interpreted as the static potential and gyM is the Yang-Mills coupling of the 
two-dimensional gauge theory. By plugging our ansatz into (I7.13P , we get: 



V 



2 L 



27c 2 a 



(#r x -i) 1 



27c 2 a 



2vrV 



An 2 a 



-(a + b). (7.27) 
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Therefore, the expected no-force condition is satisfied if a and b are related as: 



a = -b . (7.28) 

Notice that this relation solves ( 17.211) and fixes the integration constant of this equation. 
From the action fl 7 . 1 3 1) it is also possible to get the value of the two-dimensional Yang-Mills 
coupling: 

g YM [r) ATrg s 

By using the explicit expression of the function b(r) for our solution (eq. (17.251) ). and the 
radius-energy relation r = 2Tra'fi, we get: 



1 no! N c 



9Y M ( r ) 2 9s 2tt/z 2 g YM 



9ym N c 



(7.30) 



where, in the last step, we have introduced the bare Yang-Mills coupling g YM = We 
notice that, indeed, the field theory result of eq. (14.81) is reproduced by our solution. 

We finish this subsection by recalling that we have found the explicit solution for the 
functions a(r) and b(r) of our ansatz (eqs. (I7.25P and (17.281) ). while the remaining function 
H\{r, a) is given by the solution of (I7.23P (which we have not been able to solve in terms of 
elementary functions, although we think that it is doable). 

7.2 The flavored solution 

We now look for a solution corresponding to the complete fractional D1-D5 brane system, 
the inclusion of the D5-branes accounting for the presence of flavors on the gauge theory 
side. First of all, we make an ansatz of the standard form for the metric and dilaton in terms 
of two warp factors Hi and H 5 , namely: 

ds 2 = R- ZI 'E- Xl % a(i dx a dx fi + Hl /4 H 3 5 /4 (dr 2 + jXX^) + ( 7 ' 31 ) 

i=l 

e* = Hl /2 H~ 1/2 , (7.32) 
where Hi and H 5 depend on the radial variables r and a as follows: 

H 1 = Hi(r,a), H 5 = H 5 (r) . (7.33) 
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As in the unflavored case of section (17. ip . our background will contain a non- vanishing 
b = b(r) field. Actually, it is convenient to factor out in b the warp factor H$(r) and to 
represent it in terms of a new function Z = Z(r) as: 

b(r) = -Z(r)\H 5 (r)] 1 . (7.34) 

The two-forms C 2 and A 2 will have now two types of components, due to the presence of 
the flavor branes. We will represent them as: 

C 2 = (Hi 1 - l)dx° A dx 1 + C f Q 2 , A 2 = a(r)dx° Adx 1 - C f — 2 , (7.35) 

where Cf is a constant accounting for the presence of the D5-branes that wrap the orbifolded 
space and Q 2 is a two-form defined by the condition dQ 2 = |co>i A u 2 A oj 3 . The constant Cf 
is proportional to the number of flavors Nf. Actually, by imposing the quantization of the 
RR field strength F% we can write it as: 

Cf = -2g s a'N f . (7.36) 

The value of the three-form field strength T 3 for this flavored ansatz is: 



a' + (ZH^)'{H^ -I) -C f ^ * u } (dr A dx° A dx 1 + T g X 5 u l A u 2 A oj 3 ) 



r 3 H l H 5 



r 6 HiH 5 , 

(7.37) 

From the action (I7.4p it is easy to find the equations of motion. After including the Dl- 
and D5- brane sources (determined by the boundary actions (I7.13P and ( I7.14p ). we get the 
following two equations involving the NS-NS twisted scalar b: 



a' = -b' + 



Z C f 



r 3 Hl H l {r,a = 0) 



r 



3 2vr 2 a' V 4 



which are the flavored generalization of (17.211) and (17.221) . Using this result, we get from the 
equation of F 3 : 

d r (r 3 a 3 H[)' + ^(rVifgFi) = -b ,2 H 5 r 3 5 A - A^^s • (7.39) 
For any function f(r, a), let us define the Laplacian operator A as: 



A/(r,a) ^ \d r (f'r 3 ) + ^d a (fa 3 ) . (7.40) 



Then, we can rewrite (17.391) as: 

rVAffj = -b' 2 H 5 r 3 5 4 - N c K 2 orb T x 5 s . (7.41) 
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Using the relations fl7.38p - fl7.4ip one can demonstrate also that the dilaton and the Einstein 
equations reduce to: 



AH 5 = -(2n) 2 g s a'N f 5 4 , (7.42) 

By combining all the equations, one can also prove that the function Zir) defined in (17.341) 
satisfies the equation: 

1 (r 3 ZJ = 32 tt 4 (a'f 9s (n c - ^ ) 5 4 . (7.44) 

The solution to eqs. (17.421) and (I7.44p is easily found to be: 

,2N C -N f 



2ttV I 1 - 2g s a' 



H 5 = 1 + g s a'-[ . (7.45) 

The twisted RR potential which is relevant for the boundary action is A 2 , appearing in 
(17. 5p . One can show that for our solution it can be taken as: 

A 2 = -bdx° A dx l - [2a'g s N f b - An 2 g s a' 2 (4iV c - N f )] tt 2 . (7.46) 

Let us now perform a probe analysis similar to the one carried out for the unflavored case. 
Accordingly, let us consider a fractional Dl-brane extended along the Minkowski directions 
and let us expand its action (I7.13P as in (I7.26p . By using the value of A 2 displayed in (17.461) 
it is possible to show that the potential V vanishes. Moreover, one can also obtain the value 
of the Yang-Mills coupling constant, namely: 

1 - Z{T) (7.47) 



9ym{ t ) 4irg s 

By substituting in the right-hand side of (I7.47P the value of the function Z(r) given in (17.451) 
one gets: 



1 ira' 2N C -N f 1 



9ym( t ) 2 9s 47T/X 2 g YM L* 2nfi 2 V' c 2 )\ 

where again r = 2iTa'jj, and the bare Yang-Mills coupling is g\u = fiy • Then, the field 
theory result of eq. (14. 8p is also reproduced by our solution when flavor branes are added. 
Again, to find the complete solution we would have to solve for the warp factor Hi in (I7.43p . 
something we will not attempt to do here. 
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8 Summary and conclusions 



In this paper we have found string duals to gauge theories in two dimensions with M = (4, 4) 
super symmetry and a large number of colors. These duals are engineered by considering D3- 
branes wrapping a two-cycle of a Calabi-Yau cone of complex dimension two. We also 
obtained the dual of the Af = (4,4) theories with flavors both in the quenched (Nf/N c — > 
0) and unquenched limit (Nf/N c finite). We checked that our solutions correctly capture 
the perturbative running of gauge couplings, despite the fact that the resulting geometries 
develop an infrared singularity. Actually, this phenomenon is typical for the gravity duals of 
field theories with such an amount of supersymmetry. 

We have also explored the holographic realization of the Higgs mechanism for our setup and 
the meson spectrum in both the Coulomb and Higgs branches of the theory. As an alternative 
string dual we have found supergravity solutions representing a system of fractional D1/D5- 
branes on an orbifold and we have verified that these solutions also capture the one loop 
beta functions of the gauge theory (see ref. [31] for a discussion of the relation between the 
wrapped and fractional brane approaches). 

The geometries we have found present a naked singularity in the IR. As usual, such sin- 
gularity signals some interesting infrared physics. We know that, in cases with such an 
amount of supersymmetry, the singularity can be consistently screened by an enhangon (it 
is common in this case to talk about the excision of the singularity). The enhangon is the 
locus where the sources of the background become actually tensionless and the geometry of 
the supergravity solution seems to end there: it is not possible to determine the geometry 
inside the enhangon shell. On the field theory side, this corresponds to the inability to get in- 
stanton corrections. These euclidean non-perturbative configurations contribute to physical 
correlators with a factor proportional to the exponential of minus their action: exp[— ^jjr]- 
It is easy to see that at the enhangon locus g 2 — ► oo and thus these contributions become 
important (they are of 0(1)). The fact that the gravity solution ends there corresponds to 
the fact that it does not include them. This lack of information of the gravity dual is due to 
the fact that these configurations are suppressed in the large N c 't Hooft limit, according to 
which the 't Hooft coupling N c g 2 is kept fixed. Instanton corrections are thus exponentially 
suppressed in this limit. This is reflected also in the computation of the beta function. The 
gravity solution gives the exact perturbative answer but all the non-perturbative instanton 
corrections are missing: they are suppressed in the 't Hooft limit. 

This is different from the case with less supersymmetry (four real supercharges) where the 
singularity can be actually removed and the geometry is determined up to the deep IR. The 
difference is due to the fact that in this case the relevant non-perturbative configurations 
are the so-called fractional instantons, whose action goes as ^ rj2 : their contribution is 
of 0(1) in the large N c 't Hooft limit and it can thus be captured by the gravity duals. 
In the two dimensional cases we are considering here, this would correspond to finding a 
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completely regular solution describing the dual of the (2,2) supersymmetric field theory. 
There are at least two possible scenarios to engineer such a dual using wrapped branes: D3- 
branes wrapped on a two-cycle or D5-branes wrapped on a four-cycle of a three (complex)- 
dimensional Calabi-Yau. It should be possible to explicitly determine this geometry in the 
corresponding gauged supergravity approach. The D3-brane case looks more promising for 
the de-singularization of the background: a non-abelianization of the gauge field responsible 
for the twist along the cycle where the relevant D3-brane is wrapped has to be implemented 
(this is analogous to what happens in the four dimensional case studied in [53]). We are 
working in this direction |54j . 

The three-dimensional analogue of our unflavored background was analyzed in [31]. This 
background is created by D4-branes wrapping a two-cycle of a Calabi-Yau manifold in the 
type IIA theory. By using the same techniques employed here, one can add flavor to this 
solution and check that the field theory results are matched [55] . 

Another interesting point to discuss is the validity, for our flavored system, of the super- 
gravity approximation which, as is well-known, requires the Ricci tensor to be small in string 
units. As already discussed above, the geometry is singular in the IR. Let us now explore 
its behavior outside the enhangon region. Instead of looking directly at the Ricci tensor, it 
is more useful to analyze the right-hand side of the Einstein equation (1A.18I) . Clearly, the 
energy-momentum tensor of the brane (1A.20|) gives a singular contribution localized at the 
position p = pq of the flavor brane, which corresponds to the wedge shape of the function z 
(see figure [2]) . On the contrary, in the RR five-form the flavors contribute with a Heaviside 
function (see eq. (13.151) ). which is non- vanishing outside the location of the flavor source. By 
inspecting (13. 15ft one readily concludes that the contribution of Nf to F§ (and thus to the 
Ricci tensor) is just an additive constant to g. Recall that g is proportional to N c . Actually, 
one can show that in many observables N c and Nf appear in the combination iV c — which 
scales as N c in the Veneziano limit N c ,Nf — > oo with Nf/N c fixed. Therefore, outside the 
enhangon and the location of the flavor brane source, the curvature in the flavored model is 
small in the Veneziano limit, as it was in the unflavored model in the 't Hooft limit. Notice 
that this is in contrast to what happens to the flavored Klebanov-Witten model of ref. [IE] 
(see [20J for a clear discussion) and is similar to the behavior of the case studied in [12J. 

Finally, let us now comment on the validity of using the DBI+WZ action in our setup. This 
was already discussed, for example, in [21] (see section 7 of that paper). We have Nf flavor 
branes and, in principle, we would like to find a solution where they are all coincident. This 
solution would be dual to the QFT with SU(Nf) global symmetry (when we smear, we are 
breaking this global flavor group to U(l) Nf ). Notice nevertheless, that the localized solution 
mentioned above may need the DBI+WZ action to be corrected. Indeed, since the string 
coupling g s ~ N' 1 in the usual scaling for D3-branes, we find that open string diagrams, 
correcting the brane action, are weighted by factors of g s Nf ~ jf- ~ 1. This would imply 
that (unless some cancellation happens) we should not trust the tree-level DBI+WZ action. 
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Now, it is clear that the smearing avoids this potential problem. Indeed, when 'separating' 
the flavor branes, the correction to the tree-level action will be suppressed by the mass of 
the open strings between flavor branes. Of course, there is another possibility one should 
not discard, that is the fact that even in the localized case, the corrections aluded above do 
indeed cancel due to SUSY. If this is the case, the smearing should be seen as just a technical 
trick to get simpler equations. 
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A BPS equations 

We will determine the functions z(p, a), H(p, a) and g(p, a) of our ansatz by imposing that 
the background preserves eight supersymmetries. This requirement will lead us to a set of 
first-order BPS equations whose fulfillment will imply the equations of motion of type IIB 
supergravity. In order to obtain these BPS equations let us consider the supersymmetric 
variations of the dilatino A and of the gravitino ipM for the type of background that we are 
analyzing. These variations of the dilatino and gravitino are: 



5i> 



M 



2 1W ° M 

(5) F JVi-Af 5 



e* 1 



Dm + - m.N 5 T M 



(A.l) 



where $ is the dilaton, e is a doublet of Majorana-Weyl spinors of fixed ten-dimensional 
chirality and r 2 is the second Pauli matrix (which acts on the doublet e). The Killing spinors 
of the background are the e's for which 5X = 5ipM = 0. These Killing spinors generate the 
supersymmetries preserved by the background. Their existence imposes some non-trivial 
constraints on the different fields of type IIB supergravity. For example, it is clear from 
the first equation in ( 1A.1I) that the dilatino equation 5X = requires that the dilaton $ 
should be constant. Accordingly, in what follows we will take $ = 0. Moreover, the different 
components of the gravitino equation 5ipM = will lead us to the desired system of BPS 
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equations when a set of projections are imposed on the spinors e. In order to obtain them, 
let us choose the following vielbein basis for the metric (12. lh : 



e ' 1 = H-Ux ' 1 , e 2 = H- l ^d9 , e 3 = H~z ^ sin0, 

m m 

e 4 = —= da , e 5 = — ^ a (dip + cos 9 d<p) , e 6 = if 3 dp , 



e 7 = pH^doti , e 8 = pif 4 sinai<ia2 , e 9 = pif 4 sinai sina 2 ^3 • (A. 2) 

The ansatz for F5 in this unflavored case has been written in (I2.6P in terms of the magnetic 
component JF 5 = <iC 4 and of its Hodge dual. From the expression of C 4 in (12. 7p we can readily 
compute the value of JF 5 . In terms of the vielbein one- forms (1A.2I) the five-form JF 5 is given 
by: 



Th = llyl (g^e 4 + g'e 6 ) Ae 7 Ae 8 Ae 9 Ae 5 + ^ 9 f - e 2 A e 3 A e 7 A e 8 A e 9 , (A.3) 
p d <7 V / p^z 

while its Hodge dual can be written as: 



f 5 = (-g^e 6 +g'e 4 )Ae°Ae 1 Ae 2 Ae s + m29 ! 1 % ° Ae 1 Ae 4 Ae 5 Ae 6 . (A.4) 

p d <7 V / p d Z 



We can now evaluate the right-hand side of the gravitino variation in flA.lj) . In order to do 
so, we shall impose to the spinor e the following set of projections: 

T 2 345 e = e , r i23 (iT 2 ) e = e , (A. 5) 

where TmiMz- denotes the antisymmetrized product of constant Dirac matrices and the 
indices in (1A.5I) correspond to the frame (1A.2|) . After imposing these projections, one can 
show that the different components of the gravitino equation 5ipM = are satisfied if the 
system (I2.9I) of first-order BPS equations for z, H and g holds. Moreover, when the equations 
in (1 2. 9 1) are satisfied one can verify that one has the following expression for the Killing 
spinors: 

e = H~s e ~ 2 r23 e— rer e~ r78 e~ Va9 r] , (A.6) 
where rj is a constant spinor that satisfies the same projections as in (1A.5|) . namely: 

T 2 345 ^ = V » r i23(«r 2 )r? = T) . (A. 7) 

Notice that, after imposing these conditions to 77, we are left with eight preserved supersym- 
metries, as it should for a (4,4) two-dimensional gauge theory. 
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Notice that, according to its Bianchi identity, F 5 should be a closed five- form, namely one 
should have: 

dF b = . (A.8) 

The fulfillment of ( 1A.8|) is not automatic for our ansatz. Indeed, one can check that ( 1A.8I) is 
equivalent to the following partial differential equation (PDE): 



zg' 



L ap 3 H 2 J 



+ d a 



z 2 g 

ap 3 H 2 J 



m A g a 

p 3 z 2 



(A.9) 



One can verify that the equation (I A. 91) is satisfied as a consequence of the system (12.91) . Let 
us check this fact by explicit calculation. First of all, by using the third and fourth equations 
in the system (12.91) we can rewrite the two terms of the left-hand side of ( 1A.9I) as: 
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Then, by explicit calculation one can easily prove that: 
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(A.10) 



(A.ll) 



After using the first and second equation in f]2.9f) . we can rewrite the previous equation as: 

n 3 



r z g' ~ 
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ap" 



H 1 



(A.12) 



and ( 1A.9j) can be proved by using again the fourth equation in (12.91) . Thus, F 5 should be 
represented, at least locally, as dC^. A possible value of this four-form potential has been 
written in eq. (12.121) . 

Finally, one can also verify that the Einstein equations are satisfied as a consequence of 
the first-order system (I2.9P (see subsection IA.14 where this fact is checked for the backreacted 
metric). 



A.l Flavored BPS equations 

For the backreacted background, our ansatz for the metric is just the same as in the unfla- 
vored case (eq. (12. ip ). Moreover, the five-form contains an extra term (see eq. (I3.15P ). 
which takes care of the fact that the Bianchi identity is violated as in (13. 14[) . The actual 
form of T<o in the present case is: 

JF 5 = [ g da + g'dp] A u 3 A (#+ cos 9d(f>) + [g- 2vr g s (a'f N f Q(p- Pq)]v 3 A uo 2 ■ (A.13) 
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By computing the Hodge dual of JF 5 , we get: 

= 5 Z \ TT0 {g'da - zg dp) A dx° A dx 1 A u 2 + 
m l ap A ri l 

2 

(Tin 

+ -g-j [# - 2tt # s (a') 2 N f Q(p - p Q ) ] dp A da A dx° A ds 1 A (t^ + cos 8d(f>) . (A.14) 

We can now repeat the analysis performed in the unflavored case to arrive at the system 
(12. 9ft . We will still impose the projection conditions ( 1A.5I) . Actually, it is clear by comparing 
( 1A.13I) and (I A. 141) with their unflavored counterparts (1A.3I) and (1A.4I) that the resulting BPS 
equations in this backreacted case can be obtained from (I2.9P by substituting: 

g - g — 2tc g s (a') 2 Nf Q(p — p Q ) , (A.15) 

in the terms in which g is not differentiated. This is, indeed, the difference between the 
flavored and unflavored BPS systems written in eqs. (12. 9p and (13.161) . respectively. Moreover, 
the Killing spinors are still given by (IA.6p and (1A.7I) and, thus, the number of unbroken 
supersymmetries is also eight for the backreacted backgrounds. 

The modified Bianchi identity (13.141) for this flavored case is satisfied if the following PDE 
for the functions of the ansatz holds: 



9 P 



zg' 



+ d a 



. ap 3 H 2 J L ap 3 H 



z 2 g 



am 

p 3 z 2 



4 -[g -2ng s (a') 2 N f Q(p-p Q )] . (A.16) 



Eq. (1A.16P can be verified by direct calculation, following the same strategy used to prove 
(IA.9P from the BPS system. However it is simpler to check that, after subtracting an 
appropriate piece to F$, the result can be represented in terms of a four-form potential. 
Indeed, by using the equations in (I3.16P one can verify that: 

F 5 + 2ng s (a') 2 N f Q(p - p Q ) u 2 Au 3 = dC 4 , (A. 17) 

where F 5 is given by the ansatz flA.13P ~ flA.14p and C4 has the same expression as the one 
written in ( 12. 12ft for the unflavored case. It is now obvious from f 1A.17[) that F5 satisfies 

(JSH]). 

The Einstein equations in the Einstein frame are: 



Rmn — -^GmnR — — Fmpi-p* F~n 1 4 + Tmn , (A. 18) 

where Tmn is the energy-momentum tensor of the smeared flavor brane, defined as: 



Tmn -~7TuJg^- (A - 19) 
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Taking into account the form (13.131) of Sdbi we get the following expression for T MN in flat 
components: 



T 



MN 



(2tt) 4 ^(«') 2 r„_ _ I o| 1 1 Q-- -O-- 

VMN I 1 ' I — "T^T ll MP!-P 5 il NR!-R 5 V V 



(A.20) 



where f2 is the smearing form defined in (13.71) and | f2 1 is written in (I3.12p . In order to 
evaluate the different components of the tensor (1A.20I) it is convenient to write Q in terms 
of the frame forms. We get: 



n 



N f m 2 



8tt 3 Z p 3 VH 

From this expression we obtain the modulus of Q, namely: 



e 2 A e 3 A e 6 A e 7 A e 8 A e 9 



(A.21) 



N f m 2 



W \z\p 3 VH 



8(P- Pq) ■ 



(A.22) 



Then, the explicit values of the different components of T^n are: 



-Too = Tu = T 44 = T 55 = 7r g s (a ) Nf 
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\z\p 3 VH 



Kp~ pq) 
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(A.23) 



By using these values, it is straightforward (but tedious) to verify that the Einstein equations 
( IA.18I) are satisfied as a consequence of the first-order equations (13.161) . 



B The dual of the (4,4) theory from gauged supergrav- 
ity 

Let us consider the five-dimensional gauged supergravity of ref. [25] which, apart from the 
five-dimensional metric g$, contains three U(l) gauge fields A 1 , (i = 1, 2, 3) and two scalars 
ipi and y?2, which we arrange as a two-component vector (p = ((pi,<p 2 )- We also consider 
three two-component vectors a\, a 2 and CI3 whose inner products are given by: 

4 

Si ■ tij = 45ij - - . (B.l) 
According to [22] we will represent these vectors as: 

ai = (^=,V2), a 2 = (4p-V2) , a 3 = (--^,0). (B.2) 
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Then, the bosonic part of the lagrangian density of the system is given by: 
C = ^fetg~ 5 [R ~ \ {d^if - \ {d^ 2 f + Ag 2 £ 

i 

_ _ ST^ S v ip ( pi \ 2 , _ uupa\ pi p2 A 3 



(B.3) 



where i? is the Ricci scalar for the five- dimensional metric g§, F % are the field strengths of 
the abelian gauge fields A 1 ^ and g is a coupling constant. Following [251 [22], let us introduce 
the quantities X % and Xj as: 



X 1 



e 2 



- ea 



(B.4) 



This theory contains two dilatinos Aj ii = 1, 2) and a gravitino ?/> M , whose supersymmetric 
variations are given by: 



3ig 



3z 



where the quantity Vi is defined as: 



Vi = -(1,1,1) • 



(B-5) 



(B.6) 



B.l Background with (4,4) SUSY 

Five-dimensional gauged supergravity is the right theory to describe the supergravity solution 
corresponding to wrapped D3-branes. As first pointed out in [22J, the D3-brane can be 
regarded as a domain wall of the five- dimensional space and supersymmetry is realized 
by switching on appropriate gauge fields, which implement the corresponding topological 
twisting. Accordingly, let us adopt the following ansatz for the five- dimensional metric: 



21 1 



ds 2 5 = e 2 f [dxl, + dr 2 ]+— 2 [(d6) 2 + sin 2 9 (d<t>) 2 } 



(B.7) 



with / = f(r) and h = h(r) being functions of the radial variable r to be determined. We 
will consider a truncated version of the theory just described, in which the scalar field <pi = 0. 
Moreover, we will redefine the remaining scalar field tpi as follows: 



<P 



(B.8) 
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and we will assume that (p only depends on the radial coordinate r. Notice that in this 
truncated version of the theory the quantities X % and Xi defined in ( IB. 4 1) reduce to: 



e , e 



-2^ 



X 1 



In addition we will assume the following values of the U(l) gauge fields: 



A 1 = 



A 







A s = — cos 
m 



(B.9) 



(B.10) 



where m is a constant and, for simplicity, we will fix the coupling constant g as g = m. 

We shall require that the solution given by our ansatz is super symmetric, which is equiv- 
alent to demanding that SXi = 5ip^ = 0, i.e. that the right-hand side of (1B.5j) vanishes for 
some Killing spinors e. Actually, we shall impose the following projections on e: 
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After plugging our ansatz (IB. 71) . ( IB. 81) and ( IB. 101) for the metric and fields into the supersym- 
metry variations (IB.5I) . and by imposing the projections (IB. Ill) , we arrive at the following 
system of first-order BPS equations [22] : 
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2 [e 2ip - e-*) J , (B.12) 

where the prime denotes derivatives with respect to the radial variable r. In order to inte- 
grate this system, let us notice that we can get from it the following equations satisfied by 
combinations of the functions of the ansatz: 



2h' + if' = me 
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(B.13) 



Moreover, let us define a new radial variable r, related to r by means of the equation: 
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;b.i4) 



and the two new functions Ai and A 2 as: 

A 1 = 2h + V , A 2 = 2f + p- (B.15) 

Then, one can easily verify from (IB.12j) that the differential equations for Ai^ are: 



dAi 

~dT 
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dr 



(B.16) 
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whose integral is just: 



e Al « = ae 2T - I , e Aa(r) = (3 e 2r , (B.17) 



with a and /3 being constants of integration. Moreover, the equation for (p can be written in 
terms of Aj as follows: 

^ + V = I e - Al + 2 (R18) 

ar 3 3 3 

which can also be easily integrated. In terms of the original functions /, ft, and <p of our 
ansatz, the solution of the BPS system is thus: 

e 2h +V = a( ,2r _ 1 ^ 



a e 2T — - 



2 

where 7 is a new constant of integration. 



B.2 The uplifted metric 



By using the equations written in [25] one can get the expression of the metric and five-form 
of the ten-dimensional background corresponding to the solution of five- dimensional gauged 
supergravity just found. Actually, the uplifting formula for the metric is [25] : 



r/s 2 



A dsl 
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-= J2 X * ( d tf + + 9 A l f) , (B.20) 

* 8=1 



where g is the same coupling constant as in flB.3j) and, as before, we will take g = m. In 
(IB. 201) the 4> l are three angles varying between and 2n and the quantities ^ satisfy the 
relation: 



(B.21) 
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while Xi are defined in terms of the scalar field as in ( IB. 91) . The quantity A appearing in 
(lB~20l) is defined as: 



A 



8=1 



Let us parameterize the /Xj's in terms of two angles (0,ip), < 6 < tt/2, as: 

/ii = cos 6 smip , yU 2 = cos 6 cos ^ , ^3 = sin 9 



(B.22) 



(B.23) 
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When only a scalar field is non-zero as in (IB.9j) . A reduces to the following expression: 

A = e'^ cos 2 9 + e 2lp sin 2 9 . (B.24) 
Taking into account that the following combination: 

dnj = (dip) 2 + sin 2 ^^ 1 ) 2 + cos 2 ^(# 2 ) 2 , (B.25) 
represents the line element of a three sphere, one can check that : 

W + rfdtf) + e- 2 ^{dfi 3 ) 2 = e~ v A(d9) 2 + cos 2 9 dtt 2 3 . (B.26) 



8=1,2 



Taking this into account, the uplifted ten-dimensional metric ( IB. 201) when the gauge fields 
are taken as in (IB.lOj) . can be written as: 
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e~ v A(d9) 2 + cos 2 9 dil 2 3 + sin 2 9e- 2 ^ (d^ 3 + cos 9 dcj))' 

(B.27) 

Let us rewrite this metric by using the explicit form (IB. 71) of the five-dimensional metric. If 
we denote tp = (f) 3 , we obtain: 
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e'f A(d9) 2 + sm 2 9 e~ 2lp {di) + cos 9d<pf + e v cos 2 9 dVt\ . (B.28) 



B.3 The five-form 

By using the formulae written in [25] we can also get the explicit expression of the RR 
five-form F 5 . It turns out that F 5 can be written in terms of a magnetic part JF 5 and its 
Hodge dual, and that the JF 5 component can be obtained from a potential four-form C4, 
whose explicit expression takes the form (12. 7p with: 



_ e~ v cos 4 9 
9 ~ m 4 A 

Actually, the explicit expression of JF 5 is given by: 



(B.29) 
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where u 3 and uo 2 are the same as in (12.81) and (12.131) . The quantization condition of the RR 
five-form F5 is: 

F 5 = N C T 3 . (B.31) 



2 K io JM 5 
Taking into account that: 

2k\ q = (2vr) 7 ^ (a') 4 , T 3 = - 1 , (B.32) 

(2tt) 3 # s (a') 2 

and thus: 

24, T 3 = (27r) 4 g s (a') 2 , (B.33) 

we can convert (IB. 311) into: 

/ F 5 = (27r) 4 g s (a') 2 N c . (B.34) 
Jm 5 

We can use the form of JF 5 written in (IB. 301) in the flux quantization condition (IB. 311) to 
determine the value of m. To do this we must integrate JF 5 at r — > 00 along the transverse 
five-sphere parameterized by 9, 03 and the transverse three-sphere whose volume element is 
U3. It is clear that the only contribution comes from the first term in (IB.30j) . Taking into 
account that tp — > as r — > 00, one has that A — > 1 and, thus: 
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JF 5 , = — sin^ cos 3 9 d9 A d<fr 3 Aw 3 , (B.35) 
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and, therefore: 



4 (f . . „ ^ r 4v]-3 
jr 5 = / sing cos 3 / dfo u 3 = — . (B.36) 

By using this value in the quantization condition (IB. 311) . one gets that m is given by: 
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. - 47r(a') 2 g s N c . (B.37) 



Notice that m is the same as in (12.21) and that m 4 is just equal to what we called go in the 
ten-dimensional approach (see eq. (I2.20p ). 



B.4 Identification with the lOd variables 

Let us now try to identify the metric (1B.28|) with the one in (12.11) obtained in the ten- 
dimensional approach of section [2j By identifying the parts corresponding to the Minkowski 
space and the cycle one gets: 

= VA e 2f , ff-s z = VK e 2h , (B.38) 
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which leads to the following identification of the function z: 



(B.39) 



Moreover, by comparing the term of the fiber ip and that of the three-sphere dQ 2 in ( IB. 281) 
and (12. II) . we get: 
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By using the expressions of H and z obtained in (1B.38[) and (1B.39[) in (1B.40[) . we can obtain 
the expression of the two radial coordinates a and p of the ten-dimensional approach in 
terms of the variables of gauged sugra, namely: 
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In order to completely identify the two metrics (1B.28j) and (12.11) . let us calculate dp and da 
from (IB. 411) . One gets: 
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Moreover, from the the BPS system f!B.12j) , one can easily demonstrate that: 
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Plugging this result into ()B.42j) we obtain: 
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da = e f+h+v sin 9 dr + cos OS . 

m 



[BAA) 



It is now easy to show that: 



H* (dp) 2 + — (da) 2 

z 



e~ v (d9) 2 + VK e 2f (drf 



;B.45) 



from which it follows that (IB.28I) and (12.1 p are equivalent. 
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B.5 The solution in the lOd variables 

Let us now rewrite the solution (IB. 191) found by integrating the first-order system of gauged 
supergravity equations in terms of the variables used in the ten-dimensional approach. We 
will verify that the result is just that written in eqs. (12.181) . (I2.2ip and (12.221) . Actually, 
from the relation (IB. 391) between z and the metric functions h and /, it is simple to relate 
t and z. One gets: 

e 2T = 9 , 1 7n • (B.46) 
2(a — (3z) 

By using this result one can easily write the right-hand side of (IB. 191) in terms of z. It is 
convenient to write this result in terms of a new constant z*, defined as: 



One obtains: 



2 z* 



(B.47) 



-3ip 



2 z* 



z* + (2* - z) [ log(z* — z) + k] 



'B.48) 



where k is the following combination of the integration constants (3 and 7: 

k = log(2/3) - 2 7 . 



;b.49) 



By using the solution (1B.48|) in (1B.41[) . one can also get the values of cos 9 and sin# in terms 
of p, a and z. One obtains: 



cos 2 9 = 2m 2 (z* — z) p 2 



sin 2 9 



2m 2 (z t - z) 
Viz) 



(B.50) 



where the function T(z) has been defined in (12.171) . By using now the fact that cos 2 #+sin 2 9 = 
1, one can immediately obtain the equation that implicitly determines z as a function of 
(cr,p), namely (I2.18p . It is also easy to get A in the new variables: 



A 



2m (z* - z) 



P 2 + 



r 2 ( 



a 



T(z) 



(B.51) 



Using this result one can obtain the warp factor H from the identification (IB. 381) and, 
similarly, one can get the five- form function g from (IB. 291) . The result is just the one written 
in eqs. ( |Z2I) and (12T22D of section fj 
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Notice that, by differentiating the implicit relation (12.181) . one can obtain the partial 
derivatives of z(p, a). For the first derivatives one gets: 



2p(z* - z) 



p 2 + 



2a (2* - z) 



ra(«) 



a 



2 • 



p 2 + 



T(z) 



(B.52) 



One can also obtain the second partial derivatives of z(p, a) by differentiating (1B.52j) . Using 
this result one can verify that, indeed, the function z(p, a) defined in (12. 18ft solves the PDE 
written in (I2.1ip . Moreover, by plugging the values of z' and z written in (IB. 521) in the first 
and second equations of the system (I2.9P one can obtain the functions g(p, a) and the warp 
factor H(p,a). The result is, again, the one written in eqs. ( 12.211) and ( 12.221) . As a final 
check of these expressions one can verify by computing the derivatives of H and g (using 
again (IB. 521) ) that the last two equations in the system (12. 9p are also satisfied. Therefore, one 
concludes that (I2.18p . (12.211) and ( 12.221) provide the sought after solution of the first-order 
BPS system ( 1211 . 



C SUSY embeddings 

The embeddings of the D3-brane that preserve the supersymmetry of the unflavored back- 
ground are those which satisfy the condition 



r K e = e, (C.l) 

where T K is a matrix which depends on the embedding of the D3-brane and e is a Killing 
spinor of the background. In order to specify the precise form of T K , let us define the induced 
Dirac matrices on the D3-brane worldvolume as: 

7a = d a x M eS rv f , (c.2) 

where E^ denotes the vielbein coefficients, which are the ones needed to express the frame 
one-forms e M of the ten-dimensional geometry, written explicitly in ( 1A.2I) . in terms of the 
differentials of the coordinates, namely: 

e M = E™dX M . (C.3) 

Then, when the worldvolume gauge field F is zero, the matrix T K for the D3-brane takes the 
form [57] : 

r K = I -7=L=f= e a -^ 7ai ... a4 (zr 2 ) , (C.4) 
41 V-detG 4 

where 7 ai ... a4 denotes the antisymmetrized product of the induced matrices ( 1C2[) and G 4 is 
the induced metric on the D3-brane worldvolume. Let us now assume that we choose, as in 
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(13. 2p . x°, x 1 , ^ and er as worldvolume coordinates. The kappa symmetry matrix (10.4[) for 
this choice of coordinates becomes: 

1 



a/— det G 4 



(C.5) 



We will now restrict ourselves to embeddings which are of the form (I4.29p . i.e. in which the 
only non-trivial scalars are 9 and 0, which could depend on the other two CY 2 coordinates ip 
and a. The induced gamma matrices for such embeddings can be obtained from ( 1C.2I) and 
are given by: 

73.0,1 = n~\ r ,i , 



H 



li> 



la 



rn 



d^6T 2 + smOd^T 



cos 6 d,, 



tys 



<9 CT #r 2 + sin 0^1^ 



+ 



Z L 



T 4 + a cos 6* <9 CT r £ 



(C.6) 



In order to find the embeddings of the type (14.291) that are kappa symmetric and preserve 
the same supersymmetries as the background, we should compute the action of the anti- 
symmetrized product 7 x o x i^ CT (zr 2 ) on the Killing spinors e (see eq. ( 1C.5I) ). To perform this 
calculation it is interesting to realize that, from the projections satisfied by e (eq. (1A.5|) ). 
one has: 

Toi23 (ir 2 ) e = -r i45 (ir 2 ) e = e , 



r i25 (it 2 ) e = -r i34 (ir 2 ) e = -r 35 e 



r i35 ("2) e = r 0124 (ir 2 ) e = r 25 e • 
By using these conditions, one can easily verify that: 

^ x Q x ^ a {iT 2 )t = [cj + C3 5 r 3 5 + C2 5 r 2 5] 6 , 

where the c coefficients appearing on the right-hand side of (1C.8j) are given by: 



Cl 



C35 



C25 



ex H 2 £ 

[l + cos^c^] H — sin# 

^ ~ m 2 H2 



sin 9 

m 

1 

m 



a 



d^d) H 



<9 CT # + cos# (d-tpcpda-O — d^pO d a (fi) 



d^9 — a sin 6 d a c 



(C.7) 
(C.8) 



(C.9) 
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To realize the kappa symmetry condition without imposing any extra projection on the 
spinor e, the right-hand side of ( 10.8ft should contain only the term proportional to the unit 
matrix. Therefore, we must impose the following BPS conditions: 



c 35 — c 25 — , 

which reduce to the following system of equations: 



(CIO) 



sin 9 d. 



a 



d a 9 + cos 9 ( d^4> d a 9 — d^9 d a <p ) 



. 



d^9 - a sin9d a = . (C.ll) 



Notice that, if 9 and are constant, the system (IC.lip is automatically solved. This is 



precisely the embedding of flavor branes in the Coulomb branch used in section |3] (see eq. 
(13. 3ft ). Moreover, one can check that, if the BPS system (IC.lip holds, one has: 



-detG 4|BP5 = H-2 CllBPS , (C.12) 

and, indeed, the kappa symmetry condition T K e = e is fulfilled when the embedding (14.291) 
solves the BPS system ( IC.lip . Before finding its more general solution, let us try to solve 
this system by means of a more restrictive ansatz of the form: 

9 = 9(a) , = 00/0 . (C.13) 



Notice that the second equation of the BPS system (10.111) is automatically satisfied by our 
ansatz (10.13)) . while the first equation reduces to: 



sin 6* + a cos 9 o a 9 

Consistency of the above equation with the assumed dependences of 9 and in ( 10.13)) implies 
that both sides of the equation must be independent of both a and ip. Accordingly, let us 
write: 

d i ,(f> = p, (C.15) 
with p being constant. This equation can be integrated as: 

= P^ + 0o , (C16) 
where 0o is a new constant. Moreover, the equation for 9(a) becomes: 

a(l + pcos9) 
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which can be straightforwardly integrated, namely: 



° = ~ 777^ ri r , (C18) 





l+i 

P 




1 


_ sin (|) ' 




_ cos (f) ^ 





where C is another constant of integration. Interestingly, the result we have found in ( 1C.16I) 
and (1C.18I) can be compactly written in terms of the two complex variables £i and (2 defined 
in ( 14.301) . In fact one can prove that these two equations are equivalent to the following 
complex equation: 

CiCf = constant . (C.19) 

Actually, it can be proved easily that any holomorphic relation between (j and (2 solves the 
system ( IC.llj) . Indeed, let us assume that £1 and (2 are related by: 

Ci = m) , (C.20) 



where /(C2) is an arbitrary holomorphic function of £2- In order to check that fjC.20j) solves 



fjCHl]) let us take, as in (14.291) . the coordinates i[) and a as independent variables to pa- 
rameterize the embedding. Then, by computing the derivatives with respect to them of the 
relation (1C.20j) . we get: 



d* Ci = /'(C2) C2 , d a Ci = /'(C2) d„ (2 ■ (C.21) 

By eliminating /XC2) from the above equations, we arrive at the following relation: 

d$ logCi da logC2 = d^p log(2 d a logCi , (C.22) 

which, after using the definition of £1 and (2 in (I4.30p . can be shown to be equivalent to the 
system (IC.llj) . It is now clear that flC.20j) is the generalization of the solution flC.191) . 
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